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Abstract 

Critical dynamics in film geometry is analyzed within the field-theoretical approach. 
In particular we consider the case of purely relaxational dynamics (Model A) and Dirich- 
let boundary conditions, corresponding to the so-called ordinary surface universality class 
on both confining boundaries. The general scaling properties for the linear response and 
correlation functions and for dynamic Casimir forces are discussed. Within the Gaussian 
approximation we determine the analytic expressions for the associated universal scaling 
functions and study quantitatively in detail their qualitative features as well as their various 
limiting behaviors close to the bulk critical point. In addition we consider the effects of 
time-dependent fields on the fluctuation-induced dynamic Casimir force and determine ana- 
lytically the corresponding universal scaling functions and their asymptotic behaviors for two 
specific instances of instantaneous perturbations. The universal aspects of nonlinear relax- 
ation from an initially ordered state are also discussed emphasizing the different crossovers 
that occur during this evolution. The model considered is relevant to the critical dynamics 
of actual uniaxial ferromagnetic films with symmetry-preserving conditions at the confining 
surfaces and for Monte Carlo simulations of spin system with Glauber dynamics and free 
boundary conditions. 



Dynamic critical phenomena - confined geometry - finite-size scaling - magnetic properties 
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1 Introduction 



The microscopic understanding of collective dynamic phenomena in condensed matter poses one 
of the most difficult challenges for statistical physics. Accordingly the theory of these phenomena 
is in a significantly less mature state than for static properties in thermal equilibrium; also 
the corresponding experimental knowledge is very limited. At present theoretical insight into 
collective dynamics can be gained either by simulations or by studying numerically or analytically 
rather simplified models of actual condensed matter systems. Whereas in the former case the 
limitations on system size and time scales are very severe, in the latter case one has to be careful 
in accounting within the model for all those aspects of the actual systems which are relevant for 
the collective behavior under study. Understanding the link between the microscopic physical 
parameters of the system and those defining the effective dynamic models is also a crucial issue. 
Due to these difficulties only in few cases one is able to provide theoretical predictions that can 
be quantitatively compared with experiments - provided that those can be carried out in the 
first place. Nevertheless there are instances in which a universal collective behavior emerges 
which is largely independent of the microscopic details of the system and, as a consequence, also 
of the specific model used to describe it. These highly valuable circumstances arise naturally 
upon approaching a critical point, where the system undergoes a continuous, i.e., second-order 
phase transition. For universal critical properties such as critical exponents, scaling functions, 
and amplitude ratios, one is usually able to provide theoretical predictions that can be tested 
quantitatively by comparison with experimental data. In view of the universality of the critical 
properties, which is justified by the framework of renormalization-group theory and supported 
by experimental evidence, it is possible to study the collective behavior in terms of suitable 
field-theoretical models, based on minimalistic fixed-point equations of motions following from 
Landau type fixed-point Hamiltonians. This approach has been carried out successfully during 
the last decades in order to study static and dynamic critical properties of systems both in the 
bulk and in the presence of surfaces. In many cases the agreement between such field-theoretical 
predictions and (mainly Monte Carlo) simulations or experimental data is striking. 

The collective behavior of a system close to its critical point can be described in terms of the 
order parameter whose actual nature depends specifically on the system. Indeed, as long as one 
is interested in its behavior at length and time scales much larger than the microscopic ones, 
an effective Hamiltonian can be used which reflects the internal symmetries of the underlying 
microscopic system and which depends only on the order parameter and potentially a few other 
slow modes. 

Within this framework one can determine the actually observed non-analytic behavior of 
thermodynamic quantities and structure factors upon approaching the critical point. Moreover 
some of the quantities characterizing such non-analyticities (e.g., critical exponents or amplitude 
ratios) turn out to be universal in the sense that they depend only on general features of the 
effective Hamiltonian such as the spatial dimension and internal symmetries but that they are 
independent of the details of the actual system. The numerical values of the universal properties 
and of the universal scaling functions characterize the so-called universality classes [1]. On time 
scales much larger than the microscopic ones it is possible to describe the dynamics close to 
critical points in terms of stochastic evolution equations for the order parameter such that 
its resulting equilibrium distribution is given by the effective Hamiltonian of the universality 
class which the system belongs to [2]. This approach allows one to compute systematically 
the non-analytic behaviors observed in dynamical quantities, e.g., in the low- frequency limit of 
the dynamic structure factor. In turn the associated universal quantities define the dynamic 
universality class. One finds that each static universality class consists of several dynamic sub- 
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universality classes which differ, e.g., by different conserved quantities, but nonetheless exhibit 
the same static universal properties. As an example, the static universality class of the phase 
transitions in uniaxial ferromagnets is the same as that of binary liquid mixtures although their 
universal dynamic behavior is captured by two different dynamic universality classes. Various 
analytical methods, in particular the renormalization-group theory, have been developed and 
applied to provide predictions for universal quantities. (See, e.g., Ref. [3] for a recent review of 
the results obtained for the most relevant static universality classes of critical phenomena in the 
bulk. A recent summary of bulk critical dynamics can be found in Ref. [4].) 

Within this framework it is possible to account for the effects of surfaces on the critical 
behavior. Indeed real systems are always bounded by surfaces or interfaces between different 
phases which break translational invariance and thus are expected to influence the physical 
properties including universal features. In particular it turns out that compared with the bulk 
the critical behavior is locally altered within a distance from the surface of the order of the bulk 
correlation length. The resulting critical behavior depends only on general properties of the 
surface and in turn it can be classified in terms of different surface universality classes branching 
from the same bulk universality class and which are in general characterized by their own surface 
critical exponents different from the corresponding bulk ones (see Refs. [5,6] for comprehensive 
reviews) . On the other hand it turns out that there are no independent dynamic critical surface 
exponents [7]. 

In addition to the local effects near surfaces, the properties of a system are influenced by 
its finite size when the correlation length £ becomes comparable with the typical sample size 
L. Depending on the specific system and its geometry this can result even in a suppression of 
the phase transition or, generally, in a shift of the critical point and of coexistence curves which 
depends on L and vanishes in the bulk limit L — ► oo [8]. The scaling behavior that is observed 
upon approaching the critical point is expected to involve L/£ and its theoretical understanding 
is based on the finite-size scaling theory [8-10]. As a consequence of confinement and boundary 
conditions fluctuation-induced effective forces on the confining surfaces arise known as thermo- 
dynamic Casimir forces (see, e.g., Refs. [10-14] and references therein). 

Thin films provide the simplest geometry for studying theoretically the effects of confine- 
ment on phase transitions; moreover they are particularly relevant experimental realizations of 
finite-size systems. Thin film are characterized by a finite width, which in the present context 
is taken to be much larger than the typical microscopic scale, and a macroscopicly large lat- 
eral extension, i.e., much larger than the correlation length. Thin films of magnetic materials, 
confined fluids, and wetting films represent specific systems with such a geometry which are 
indeed investigated experimentally. Their static critical properties have been theoretically and 
experimentally investigated in the past for different universality classes and boundary conditions 
(see Refs. [15,16] and references therein). 

Dynamics in confined geometry is, instead, a less explored subject, both at the critical point 
and below. Novel phenomena have been observed in the dynamics of phase separation [17] 
occurring in the two-phase region of the phase diagram of confined binary liquid mixtures, after 
a quench from the homogeneous state. In particular the interplay between surface-directed 
spinodal decomposition (see, e.g., Ref. [18]) and confinement has been studied numerically in 
Ref. [19] for a symmetric binary mixture with purely diffusive dynamics (Model B in the notion 
of Ref. [2]), a simplified form of the actual dynamics of fluid mixtures (Model H [2]). At 
the critical point, which is the focus of the present study, most of the theoretical results have 
been obtained for the case of a finite hypercubic geometry with periodic boundary conditions 
and purely dissipative dynamics [1, 20-25] (Model A in the notion of Ref. [2]) or dynamics 
coupled to a conserved density [26] (Model C in the notion of Ref. [2]). The dynamic structure 
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factor, the spin transport, and the thermal conductivity of the three-dimensional XY model on 
a cubic lattice with periodic boundary conditions have been studied by means of Monte Carlo 
simulations [27]. In all the cases mentioned above, due to the translational invariance, there are 
no surface, i.e., spatially varying effects. Monte Carlo simulations have been performed to study 
the thermal conductivity of the planar magnet lattice model in a bar-like geometry {H x H x L 
with L 3> H) with open boundary conditions [28], aiming for a comparison with experimental 
results for 4 He at the superfluid transition confined to an array of pores [29]. The same problem 
has been addressed within the field-theoretical approach, studying Model F dynamics [2,30] in 
aixLxoo geometry with Dirichlet boundary conditions (DBC, i.e., vanishing surface fields) 
for the order parameter [31]. In both cases the agreement with experimental data is quite 
good. Critical dynamics in the film geometry has not yet been studied systematically, in spite 
of available experimental data for some specific systems. In particular, the thermal conductivity 
of 4 He close to the normal-superfluid transition and in confined geometry has been investigated 
experimentally in some detail. Field-theoretical methods have been employed to analyze the 
so-called thermal boundary resistance (Kapitza resistance) between the superfluid 4 He and the 
wall confining the system. This can be carried out by considering Model F dynamics in a semi- 
infinite space with DBC [16,32]. In spite of these results, for a specific surface quantity, the 
theoretical prediction for the full finite-size behavior of the thermal resistance across a film is still 
lacking. Moreover, recent experimental findings [33] are in disagreement with the field-theoretical 
predictions of Ref. [32]. It has been argued that this might be a consequence of the choice of DBC 
being inappropriate for describing helium confined to a film. Some other transport properties 
have also been investigated theoretically for the film geometry. In particular, in Ref. [34] the 
effects of confinement on the critical diffusivities have been studied within the decoupled-mode 
approximation [35] for the dynamic universality class of liquid- vapor phase transitions in a one- 
component fluid (Model H in the notion of Ref. [2]) and for the superfluid transition (described 
by Model E [2]) with DBC at the confining plates. The diffusion constant considered is the 
one associated with the density current and the superfluid flow, respectively. In Ref. [36] the 
finite-size behavior of the ultrasonic attenuation that is observed upon approaching the critical 
point of the superfluid transition in 4 He has been studied. The sound velocity can be related 
to the frequency-dependent specific heat (see Ref. [36] and references therein), allowing for a 
quite direct field-theoretical analysis within the Gaussian model. Instead of the full Model F 
dynamics, which is the appropriate one to describe the superfluid transition, it is possible to deal 
approximately with this problem by considering the simpler Model A dynamics of the superfluid 
order parameter. Then, by applying Dirichlet boundary conditions, the scaling functions for 
the ultrasonic attenuation have been computed, resulting in a good agreement with available 
experimental data [36]. 

Recent efforts [37,38] address theoretically some aspects of non-equilibrium (critical) dynam- 
ics of a scalar fluctuating field 4> (which can be, e.g., the order parameter of an Ising ferromagnet 
or the deformation of an elastic membrane) in film geometry with DBC. The dynamics of 4> is 
assumed to be purely dissipative (as in Model A), the effective Hamiltonian is taken as a Gaus- 
sian, and the immobile confining walls to be actually "immersed" in the fluctuating medium. 
However, different from the usual Model A, the fluctuations of the field are taken either to be 
due to external forces [37] or as quasi-equilibrium thermal noise generated by a space- and time- 
dependent temperature profile [38]. The main focus of these analyses is the computation of the 
Casimir-like non-equilibrium fluctuation-induced force that acts between the confining walls for 
different instances of driving forces and temperature profiles, whereas little attention is paid to 
the actual dynamics of the field <j) in the space delimited by the walls. 

In spite of these results a systematic investigation of the critical dynamics in film geometry 
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is still lacking. In view of the rapidly developing experimental techniques able to resolve space- 
and time-dependent quantities on the proper mesoscopic scale it is important to provide the- 
oretical predictions for experimentally accessible quantities such, e.g, time-dependent response 
and correlation functions, dynamics of fluctuation-induced forces, etc. 

In the following we set up the field-theoretical description of the purely relaxational dynam- 
ics (Model A in the notion of Ref. [2]) in film geometry with Dirichlet boundary conditions. 
We consider the static universality class of systems with a TV-component order parameter and 
0(iV)-symmetric interactions, such as the Ising model (N = 1) or the isotropic XY (N = 2) 
and Heisenberg (N = 3) models. Although it is possible to study the relaxational dynamics for 
general N by kinetic Monte Carlo simulations, for Model A N = 1 is the only physically relevant 
case, experimentally realized in anisotropic magnets. For systems with N > 1 the actual dy- 
namics requires a description in terms of more complex models [2,39]. Moreover also the proper 
description of mixing-demixing transitions in binary liquid mixtures and liquid-vapor transi- 
tions in one-component fluids, whose static universal properties are given by the case N = 1, 
calls for different dynamical models and boundary conditions [2,17]. Accordingly, the model we 
consider here is relevant for the dynamics of actual uniaxial magnetic films (without energy con- 
servation [2]) with symmetry-preserving boundary conditions and for Monte Carlo simulations 
with Glauber dynamics of O(N) order parameter models with free boundary conditions. Keep- 
ing in mind these caveats the present analysis provides a theoretical framework which might 
nonetheless turn out to be useful also for the more complicated dynamical models. Explicit 
expressions for the universal scaling functions are obtained within mean-field theory (Gaussian 
model), whose actual behavior (independent of N) turns out to be already quite rich. We discuss 
in detail the surface behavior close to the confining walls and the temporal crossover between 
different regimes. The analysis carried out here lends itself for possible, future extensions within 
the field-theoretical approach to include the effect of fluctuations beyond mean-field theory. 

After this introduction the paper is organized as follows. In Sec. |2]we recall some general 
scaling properties which will be useful in the following. In Sec. |3] the model is described and the 
universal scaling functions for the dynamic Gaussian linear response and correlation functions 
are derived and studied in detail in Sec. HI In addition we discuss the effects of time-dependent 
fields on the fluctuation-induced Casimir force acting on the confining walls of the system, 
providing analytic expressions for the associated Gaussian scaling functions. In Sec. El the 
nonlinear relaxation of the order parameter is analyzed at the bulk critical point of the model 
and the crossover between bulk-like, surface-like, and eventual linear relaxation is emphasized. 
In Sec. we provide a summary of the main results. Most of the details of the computations 
are reported in the Appendices El and lEHUl Instead, in Appendix El the computation of the 
bulk universal amplitude ratio associated with the divergence of the relaxation time is reported, 
whereas in Appendix Owe determine and discuss the useful analytic expression for the static 
order parameter profile across the film in the low-temperature phase and for Dirichlet boundary 
conditions, which, to our knowledge, has never been reported in the literature. 

2 General scaling properties 

We consider a confined system in d dimensions with film geometry oo rf_1 x L and Dirichlet- 
Dirichlet boundary conditions (corresponding to the so-called ordinary-ordinary surface univer- 
sality class [6]) and purely dissipative relaxational dynamics (Model A of Ref. [2]). 

For future reference we introduce here some of the notations used in the following. We define 
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the reduced temperature 

T = ^ — - (!) 

J-c,b 

where T is the temperature and T c ^ is the transition temperature in the bulk. With £ and Tr we 
denote the true correlation length and the true relaxation time, respectively, as the characteristic 
length and time scales defined via the exponential decay of two-point bulk correlation functions 
in thermal equilibrium. Upon approaching the critical point (r — ► 0) both £ and Tr diverge 
with the following leading singularities: 

i{r-*Q k )=Z±\T\- v (2) 

and 

T R (r - ± ) = TtV\~" M = Ttm/tfY (3) 

where v and z are universal standard bulk critical exponents whereas £^ and are non- 
universal amplitudes depending on the microscopic details of the system. Within mean-field 
theory (MFT) corresponding to d > 4 one has v = 1/2 and z = 2, whereas, for the Ising 
universality class (N = 1) in d = 3, v = 0.6301(4) [3] and z ~ 2.02 (see Ref. [40] and references 
therein for a summary of the various estimates of z) . The values of ^ and are different for 
r — > + and r — > 0~, forming universal amplitude ratios £o~/£o~ and Tq~ /T ~ with Co"/CcT = 
within MFT, whereas £o"/£(7 = 1.896(10) 1 for the Ising universality class in d = 3. For 
the second-moment correlation length similar results can be found in the literature [3,41]. In 
Appendix the ratio /Tq~ is computed to first order in e = 4 — d for the Ising universality 
class of Model A, leading to T + /T ~ = 2 within MFT and /T Q = 3.3(4) in d = 3 2 . 
The bulk order parameter m vanishes for r — > 0~ as 

m = mof-TV 3 , (4) 

with the universal exponent (3 and the non-universal amplitude mo; within MFT (3 = 1/2, 
whereas for the Ising universality class in d = 3, (3 = 0.3265(3) [3]. 

In the following we will be concerned with quantities G defined in the film geometry. They 
generally depend on a set {x, t} of spatial coordinates and times, on the temperature (expressed 
in terms of r), and on the film thickness L. Since upon approaching the critical point the 
dominant length and time scales are given by £ and Tr, respectively, the following scaling 
behavior is expected in the critical region |r| >C 1: 

x,t};r ^ 0,L) =o± ( X) ° ({x/£, t/T R }; L/i) 



=o± l^j F^ ± ({ X /L,(t/T^/LY};L/0 

where are non-universal constants which have the same engineering dimension as the observ- 
able G and can be expressed in terms of £q~, mo, Tq~, and universal amplitude ratios. is the 

(i) 

scaling dimension of the quantity G, and i = 1,2, are universal scaling functions. The 

1 This quantity is denoted by U( g in Ref. [3]. The numerical value quoted here, obtained by combining 
high-temperature expansions with a parametric representation of the equation of state, is taken from Tab. 11 
therein. 

2 From Eq. fTTTt one has T + /r o " = 2(1 + 2e In 2/3) + 0(e 2 ). In order to gain a rough estimate of this ratio for 
e — 1 the [0, 1] and [1, 0] Pade approximants can be used, yielding the value 3.3(4). 
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second line in Eq. (J5J) is the scaling form, equivalent to the first one, in which we shall present 
our results. For the two-point correlation function C (see Sec. 01 one has Ac* = d — 2 + ij (which 
agrees with the static two-point correlation function and defines the static bulk critical exponent 
77 with rj = within MFT, whereas r\ = 0.0364(5) [3] for the three-dimensional Ising universality 
class) whereas for the response function R, Aft = Ac + z = d — 2 + rj + z. For the magnetization 
one has A m = f3/(yz) (A m = Aq/(2z) if hyperscaling holds). Crossovers between surface and 
bulk singular behaviors characterized by surface and bulk critical exponents are related to the 
singular behavior of the scaling functions if x approaches the confining walls (see Sec. EJ). 
In this respect one has to keep in mind that the scaling properties (Eq. iJSJ) only hold in the 
scaling limit, i.e., distances between two spatial points, distances from confining walls, and time 
differences must be sufficiently large compared to microscopic scales. The field-theoretical ap- 
proach has been proven to be a powerful tool to compute both the exponents and the scaling 
functions appearing in Eq. © for various measurable quantities. In the following we present 
the mean-field form (i.e., tree- level approximation in the field-theoretical language which is valid 

(2) 

for d = 4 up to logarithmic corrections) of the scaling functions for various observables. 
In many cases a reasonably good agreement between experimental or simulation data and field- 
theoretical computation is already obtained by using mean-field scaling functions combined with 
higher-order estimates for critical exponents entering into their scaling arguments. 



3 The Model 
3.1 Definition 

The time evolution of a ^-component field </?(x, t) = (v?i(x, t), i = 1, . . . , N) under purely dissi- 
pative relaxation dynamics (Model A of Ref. [2]) is described by the stochastic Langevin equation 

3^(x,i)=-Q-^Mr + C(x,i), (6) 
dp{x, tj 

where £1 is a kinetic coefficient, £( x , i) a zero-mean stochastic Gaussian noise with correlations 

(Ci(x,t)C J (x',i')) = 21M(x - jd)8{t - t%j, (7) 

and TC[p] is the static Hamiltonian. The universal properties near the critical point of a second- 
order phase transition are captured by the Landau-Ginzburg form [1,6,42,43] 



Hfr] = [dV 
Jv 



- I d^XH V, (8) 
<dV 1 



where tq cx T is a parameter that takes the value ro jC rit for T = T c ^ (ro ;Cr it = within MFT) 
and go > is the coupling constant providing stability for r < 0. The surface term implies the 
boundary conditions (p = c~~ 1 d x± p at x± = 0, L such that the fixed-point value c = 00 leads 
to Dirichlet boundary conditions as considered in the following; we do not consider surface field 
contributions h s Jgyd d_1 X| ip. We use the notations V = R^" 1 x [0, L] and d^ = d d ~ 1 x.\\dx±, 
where the position vector x = (x||,xj_) is decomposed into the d— 1-dimensional component Xm 
parallel to the confining planar walls and a one-dimensional one x± perpendicular to them. 

Instead of solving the Langevin equation for ip [Q and then averaging over the noise distribu- 
tion P[C], the equilibrium correlation and response functions can be directly obtained by means 
of a suitable field-theoretical action S[ip, <p~] [1,42,43] so that, for an observable 0[<p], 

(O) = j [dC] OM(]]P[(] = j [dpdp] Oe- s ^ . (9) 
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(Note that within the conventions we adopt, Jfdi^dc^] e s ^'^ = 1 [43].) For the Langevin 
equation © with the Gaussian noise (Eq. Q) the field-theoretical action is given by [1,42,43] 



S[(p, <p] = J tit J dV 



ifidtf + Clip ^r-^- — <pCl<p 
dp 



(10) 



where <p(x, t) is an auxiliary field, conjugate to an external bulk field h which linearly couples to 
the order parameter p so that TL[p,h] = Tt[p] — f dVhp. As a consequence, for an observable 
O the following relation for the linear response to the field h holds: 



5{0)h 



5h(x,t) 



= n(p(x,t)o) h=0 , (ii) 

h=0 



where is the average taken with respect to the action S[p, p; h] associated with TL[p, h}. In 
view of Eq. (|lljl. p(x,t) is called response field. In the following we will be mainly concerned 
with the response of the order parameter field to the external perturbation h, given by 



Sh(x.i,ti 



= fl(p(x 1 ,ti)p(x2,t 2 ))h=o- (12) 

h=0 



Causality implies that {p(x.2,t2))h does not depend on /i(xi,ti) whenever t\ > t%, i.e., the order 
parameter at a given time does not depend on possible perturbations at later times. Accordingly 
i?(xi, t\\ X2, t%) vanishes for t% < t%. 

The effect of confining walls in the case we are interested in amounts to the Dirichlet boundary 
condition for the field p (i.e., infinite surface enhancement c; see above and Ref. [6]), 

p(x B ,t) = 0, Vt, (13) 

where we denote by xg the position vector on the boundary dV. The Gaussian (i.e., go = 0) 
equation of motion for the field p> given by —dtp + Ap + r^p) = and Eq. ()13|) yield the 
boundary condition 3 . 

00^1^=0 = 0, (14) 

which is fulfilled by imposing 

p(x B ,t) = 0, Vt. (15) 

In order to diagonalize the Gaussian part of Eq. (|10j) it is useful to decompose both fields p 
and p> in terms of eigenfunctions of the Laplacian A fulfilling the boundary conditions (f 13ft and 
(fT5|) [7,11] according to 



•«)-E/ 



e i (P- x H-^)0„(p,a;)$ n (x ± ;L), (16) 



where (f> = (p,<p, and 

, s „-i(2vr) d -i 7 R 2vr V 7 

with p as the d — 1-dimensional momentum parallel to the confining walls. The transverse 
momentum takes, instead, discrete values k n = nn/L with n = 1,2,.... The eigenfunctions 
&n (x± ; L) are given by 

$n(x±;L) = y / 2/Lsm(k n x±), (18) 



3 Indeed J Q dx±tpdl J _<p = f dx±ipdl ± tp + {ipd x± <p - <pd m± tp)\x ±=0 , where f\% ±=0 = f{x± = L) - f(x± = 0). 
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so that & n (x± = 0;L) = $ n (x± = L;L) = 0. Note that for the total momentum q n = (p, k n ) 
one has |q n | > ir/L and thus the homogeneous fluctuation mode is suppressed by the boundary 
conditions. In terms of the functions introduced above, the Gaussian part So of Eq. (|TU)) can be 
written as 

So[0,<p] = ljt[ (^(-p,-^),^(-p,-^))Mf < ? n 5 P '^V (19) 

where the inverse propagator M is given by the matrix 

-iu + fi(q£ + r 

iu + n(c& + r ) -2VL 



I ■ ( 2 0) 



Within MFT, the two-point response and correlation function and C^°\ respectively, are 
determined by M : 

i? iilnin 2 (P 1 ' U;i 'P 2 ' tJ2 ) = n (^ini(Pl> w l)^2n 2 (P2,W 2 )) flo =0 

= (2vr) d ^- 1 )(p 1 +p 2 )5(u 1 +tu 2 )5 ili2 6 nin2 R^( 

where and indicate the field component and the Fourier mode according to Eq. (|16|) . 
respectively. For the correlation function one finds 



C 'i 1 i ) 2 n 1 n 2 (Pl» a; i;P2^2) = (fami (Pi, ^l)& 2 n 2 (P2, U 2 ))g =0 

i 2 &nrn 2 C (qn 2 > w 2 ) 

while (c^j ini (pi,wi)(/5j 2n2 (p2,u;2)) 90 =o = due to causality. From Eq. (|2()j) one obtains 



(22) 



R(P) ( q , w ) = _ " (23) 

-zw + ft(q 2 +r ) 

and 

90 

C(°) (q, u) = - .J" ^ . (24) 

w 2 + [Q(q 2 + r )] 2 

Note that and C<M are the mean-field response and correlation function, respectively, of 
the system in the bulk. Within MFT the presence of the boundaries is accounted for by the 
spatial dependence of the eigenfunctions & n (x±;L) (see Eq. (|T%|0 and by the quantization of 
the allowed momenta. 

For later purposes it will be useful to provide for these functions their representation in terms 
of time and transversal coordinates Xi±. Following Eq. (|16|) we define 

0(P, x ± ,t) =J2 2~e-^0 n (p, u)<S> n (x ± ; L), (25) 

71=1 ^ 7F 

and 

oo 

0(p, u) = ^2 4> n (p, u)$ n (x±; L), (26) 
n=l 

so that, according to Eq. l(2*2*|) . 

C 'i 1 i 2 (Pl ) ^l-L^i;P2,a;2±,fe) = {iPi 1 (Pl,X 1 ±,tl)(pi 2 (p 2 ,X 2 ±,t 2 ))g Q=0 ^ 

= (2TT) d - 1 6^ d - 1 \p 1 +p 2 )5 lll2 C^\p 2 ,x 1± ,x 2± ,t 2 -t 1 ) 



9 



and 



= (2n) d 5^ d - 1 \p 1 +p 2 )5(u; 1 +L0 2 )5 tli2 C^(p 2 ,x 1± ,x 2± ,u; 2 ) 

with an analogous expression for the response function. Using the above formulae we obtain 
(recalling that q n = (p, ttti/L)) 

C^(p,x 1± ,x 2± ,t) = Y,^n(xi±;m n (x 2± ;L) /_ e -^c(°)(q n ,a;) (29) 

n=l 

and 

oo 

= Y,^{xi^m n {x 2 r,L)C^\ ), (30) 

n=l 

with C(°)(q,w) given in Eq. (|2*1)) . In Appendix [B] we show how one recovers the known results 
for the equal-time correlation function in the film geometry discussed in Refs. [15] and [44] for 
various boundary conditions and in Ref. [6], Sec. IVA, in the case of a film with equal boundary 
conditions. Analogous relations hold for the response function. 

3.2 The fluctuation-dissipation theorem 

For small fluctuations within equilibrium dynamics the two-point response and correlation func- 
tions are not independent quantities. The relation between them is provided by the fluctuation- 
dissipation theorem (FDT): 

= -R(t) , for t > , (31) 

at 

where with C{t) and R(t) we indicate summarily the time-dependent two-point correlation 
function and the response function, respectively. As equilibrium quantities their dependence on 
two time variables reduces to a dependence on the time difference only. Indeed the theorem 
is a consequence of the time-translation invariance and time-reversal symmetry of equilibrium 
dynamics. Keeping in mind that correlations vanish in the long-time limit, one thus has 

/oo 
ds R(s) , for t>0. (32) 

Time-reversal symmetry in the equilibrium state implies C(t) = C(—t), which, combined with 
Eq. (|32j) . allows one to determine completely the correlation function from the response function. 
Let us recall that the causality of the response function (linear or not) implies R(t) oc 0(t) where 
0(t) = 1 for t > and otherwise. For later purposes it is useful to express the FDT also in 
other forms. Defining the Fourier transform of C(t) as 

/+oo 
dt e iujt C(t) (33) 
-oo 

and with an analogous definition for R(u>) the FDT (Eq. (|31jl) can be written also as 

C(u) = -ImRfu) (34) 
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where Im indicates the imaginary part of the expression. Equation (|32|) and causality yield 
another useful form of the theorem: 



C(t = 0) = R{uj = 0). (35) 

It is straightforward to ver ify that and given in Eqs. (|23|) and (1241) as well as 
Eqs. (|29|) and ()3(Jj) satisfies the FDT ()34j) . Moreover it can be shown that, as expected on 
general grounds, this is true also if the effect of fluctuations are taken into account, i.e., beyond 
MFT. 



3.3 Some mean-field results 

According to the general scaling properties discussed in Sec. 121 within the Gaussian approxima- 
tion one can identify the correlation length £ and the relaxation time Tr from Eqs. (|23[) and (|24|) : 

£(r > 0) = r„ 1/2 (36) 

and 

T r {t > 0) = (ftro)- 1 , (37) 

respectively. Thus within MFT v = 1/2 and z = 2 as expected, and the kinetic coefficient in 
Eq. (jHJ) can be expressed in terms of the experimentally accessible (non-universal) amplitudes 
£o and To (see Eqs. @ and (jHJ)) as 

n = ^+ (38) 

and 

-o = ^ . (39) 
Moreover, for T < T c ^ the mean-field equation of state for the bulk order parameter m leads to 

m = (-6r / 5o ) 1/2 , (40) 

i.e., f3 = 1/2, and 

90 = 6(m £ + )- 2 . (41) 
The relations in Eqs. (|36j) - (|41|) hold for the present continuum model (Eqs. ©-©)• 



4 Linear behavior 

In this section we study in some detail the behavior of the response and correlation function for 
the model introduced in Sec. |SJ By virtue of the FDT linear response and correlation functions 
are not independent but the correlation function can be obtained from the response function and 
vice versa. We shall first focus on the response function in Subsec. 14.21 and then in Subsec. 14.41 
we shall determine correlation functions by applying the FDT. 
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4.1 Scaling forms for the response and correlation functions 

For future reference we provide here the general scaling forms for some of the quantities we 
shall discuss in the following. As already stated in Sec. E] (see Eq. ©), scaling occurs upon 
approaching the critical point. In the specific case of the two-point response function in the 
(p, x±, u;)-representation one has 

R{p,x 1± ,x 2± ,u) = o%(J^J 1 k ± (pL,x l± /L,x 2± /L,LaT±(L/t±) z ,L/0, (42) 

where 6^ are non-universal amplitudes which we fix to be equal to the corresponding bulk ones. 
The functions &± are universal scaling functions. In the (p, x±, ^-representation, this reads 

R(p,x 1± ,x 2± ,t) = ^(^) " £±( P L,x 1± /L,x 2± /L,(t/T±)(^/L) z ,L/0 , (43) 

where the universal functions &± are the Fourier transforms of M± with respect to their fourth 
argument. 

Analogously, for the correlation function in the (p, x±, u)-representation one has 

C(p,x l± ,x 2± ,uj) = 6± (Jzj tf ± (pL,x 1± /L,x 2± /L,LoT±(L/t±y,L/0 , (44) 

where 6^ are non-universal amplitudes, which again we fix to be equal to the corresponding 
bulk ones, and the functions < rf± are universal scaling functions. In the (p, x±, ^-representation, 
this reads 

C(p,x 1± ,x 2± ,t) = o±4 (4) \±(pL,x 1± /L,x 2± /L,(t/T±)(t±/Ly,L/0 , (45) 

where the universal functions c £± are the Fourier transforms of c €± with respect to their fourth ar- 
gument. In view of the previous scaling forms it is convenient to introduce the following suitable 
set of dimensionless scaling variables, defined as p = pL, Xi± = Xi±/L, t^ 1 = (t/T^)(^ / L) z , 
u)± = uT±{L/£fy, and L = L/£. 

As stated above, concerning the non-universal amplitudes 6^ and 6^ we consider the corre- 
sponding correlation functions (Eqs. (|42|) and (|45|) ) in the bulk. The critical structure factor in 
the bulk is given by 

C£f(q,i = 0) = ^, (46) 

which defines the non-universal amplitude (see, e.g., Ref. [41]). Here and in the following 
using the subscript "crit" means that the function corresponds to r = 0, i.e., to bulk criticality. 
Beyond MFT (in this case two-scale universality holds [41]) can be expressed in terms of 
the universal amplitude ratios Q3, R c , and (see Ref. [41] for their definitions and numerical 
values) and the non-universal bulk amplitudes mo and ^ (see also Ref. [15]), as 



Doo = ^ml(^r^ . (47) 



From Eq. (|46|) . via a Fourier transform in one of the d dimensions, one finds 

C*T(P, xx±,x 2± = x 1± ,t = 0) = g v p- 1+r > (48) 
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where Qy has been introduced in Appendix A of Ref. [15] and is given by 



A» r(i/2 - 77/2) 
^ - rd-,/2) • (49) 



In view of Eqs. ()151) and 1)15)1 this leads to 

~± r 1 r(l/2 - r?/2) Q 3 .R C 2 ± ± 

°c = ^ r oft>) =^TF^T(4? o o(g ' (50) 



The FDT (see Eq. (|82)l) establishes the following relation between the correlation and response 
function: 

C(p, xi±, x 2 ±, t = 0) = i?(p, xi±, x 2 ±, u = 0) (51) 
so that in the bulk one has 

I&(p,xi±,x 21 = xix, w = 0) = ^p" 1 "** . (52) 

Comparing this equation with Eq. (|42)l leads to 



ft ± r ,,±a-„ 1 r(l/2-7y/2) Q 3 R C 

°* = M£o) "=2^ ra-ry/2) (i?^™ ^ } " (53) 

Thus the non-universal amplitudes of R and C are determined by the experimentally accessible 
non- universal bulk amplitudes mo, £n" , and T + . This fixes the normalization of the scaling 
functions ffl± and c €±. For the present model and within mean-field theory these non-universal 
amplitudes are 

= 1 (54) 



and 



so that 



and 



QP = \ (55) 



8*°> = (56) 

8*°> = f . (57) 

Here and in the following with the superscript (0) we indicate the mean-field value of the 
quantities which the superscript refers to. The FDT (see Eq. 1)32)1 ) provides, together with 
Eqs. (|5Uj) and (|53jl . the following relation between the scaling functions &± and c €± (and thus 
between M and ff): 

^±(p,xi±,x 2 ±,t,L) = ds^±(p,x 1± ,x 2 ±,s,L). (58) 

In the following we shall be mainly concerned with the case r > 0. In order to avoid a 
clumsy notation we thus shall omit in the following the specification ± from scaling forms and 
amplitudes. 

In Subsec. 14.21 we shall discuss the behavior of the response function in the semi-infinite 
geometry close to a confining wall (see, c.f., Eq. I)9U)I ). It can be obtained from the scaling 
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function Eq. (|43|) in the limit L — > co with £, and i fixed. Thus one expects a well-defined 
limit for the response function, i.e., 

M{pL, x 1± /L, x 2± /L, (t/T )(^/L) z , L/0 

r\-(i-v-z) _ (59) 



where M 00 / 2 is the scaling function for the semi-infinite geometry. By using the short-distance 
expansion [6] one easily concludes that for — > (i.e., <C £, |p| -1 , Co(V^o) 1 ^) 

^°° /2 (pe, x^/e, wc, (t/T )(e /o 2 ) 

where /3i is the critical exponent for the surface magnetization [6] (f3\ = 1 at the ordinary 
transition within mean-field approximation, whereas j3\ ~ 0.77(2) [6] for the three-dimensional 
Ising universality class). Considering the case p = and T — > T c ^ (i.e., £ — > co) one expects 

lim M™ /2 (0, y) = Vy-m-WivM-*!-')/* , (61) 
y— >o 



where P is a universal constant. Thus, for T = T c ^ and Xi± <C £o(*/^b) 



1/z 



Roo/2 (p = 0, XLL , x 2± , t - oo) = ^ ( ) (62) 



where is the response function in the semi-infinite geometry. 

In Sec. I4.4I we shall discuss in detail the mean-field behavior of the correlation function 
C(p,x±,x±,tu) in planes parallel to the confining walls. For this function we present in the 
following some scaling properties valid beyond the mean-field approximation. 

From a short-distance expansion one concludes that the scaling function c € (see Eq. (|44j) ) of 
the correlation function C behaves as 

V(pL, x 1± /L, x 2± /L, ujT (L/^) z ,L/0 

— ;H %w{pL,uTo{L £q) 

for <C L, ^(wTo) -1 / 2 , £, |p| _1 , where is the universal scaling function associated with 
the behavior close to the walls (in this case the one located at x± = 0). The behavior of 
allows one to define the universal constant 

A W = ^ W (0,0,0) , (64) 

that appears in the scaling behavior 

C C r it (p = 0,x ± ,x ± ,u = 0) = b c A W ( -) (65) 

of the critical (i.e., r = 0) correlation function for x± <C L. Analogously we can define the 
following universal constants: 
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lim <g w (0, v, 0) = BZv 2 ^ -WW-O-v+'V* , (67) 

v^oc 

and 

lim tf w (u, 0, 0) = c^lul 2 ^ 1 -^^-^-^ (68) 
entering into the scaling functions 

C(p = 0,a; ±) x ±) w = 0) = 6o^ f "T ( 69 ) 



for x ± < £ < L, 



C crit (p,x ± ,x ± ,u; = 0) = 6 c C^(|p|Co)- (1 -^ )+2(/3l - /3)/l ' ^ (71) 



for (£ /L) z « LoT « (£ /z±) z , and 

for 1/L <C |p| <C respectively. 

From the previous equations we recover the values of well-known surface critical exponents 

(s) (s) (s) 

for the semi-infinite geometry, i.e., ov , cr^ , and <jp [7,45]. They describe the divergence of 
the two-point correlation function parallel to the surface, so that for p = and to = it diverges 
~ r °" T , for r = and u = it diverges ~ |p| CTp , whereas for p = and r = it diverges for 

(s) 

to —> as uj~ a ^ [45]. 

As far as the behavior of the correlation function in the middle of the film, i.e., C(p, L/2, L/2, ui) 
is concerned we define (see Eq. (|44|)1 

<r( P L,l/2,l/2,u,T (L/£ ) 2 , I/O = ^/(pL, uT (L/ £ ) 2 , L/ £) (72) 

where ^/ is a universal scaling function. As already discussed for we can define also from 
^/ the following universal constants: 

%(0,0,0) = A 1 , (73) 

lim SfAO, 0, w) = ALw^ 1 -^ , (74) 
w— >oo 

lim tfj(0, u, 0) = Bl.v-Q-- 1 *-*)/* , (75) 

v— >oo 

and 

lim ^/(u,0,0) = C4 > |ur (1_T?+2) • (76) 

u— >oo 

These constants enter into the following scaling functions: 

1 — Tj+Z 



C crit (p = 0, L/2, L/2, cu = 0) = -4 7 o c f |^ 

C(p = 0, L/2, L/2, lo = 0) = o c ^ (1) " + for £ « I 



(77) 
(78) 

C crit (p = 0, L/2, L/2, w ) = ooB^^To)-^-^)/ 2 for u,T « (|) , (79) 
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and 

C CTit (p,L/2,L/2,u; = 0) = bcC I OD (\p\^r {1 ' v+z) for |p| > j , (80) 

respectively. In Sec. 14.41 we shall confirm these scaling forms within mean-field approximation 
and determine also the mean-field values of the universal constants involved. 

4.2 Response Function 

Our aim here is to discuss, in different representations, the response function introduced in 
Sec. 01 Combining the analogue of Eq. (|29[) for the response function and taking into account 
the explicit expression in Eq. (|23|) we have 

/+OOJ 
^ e --'i?(°)(q n ,u,) = 9(t)ne- n ^ +r ^ (81) 
-oo 27r 

and thus 

oo 

R {0) (P, x 1± ,x 2± ,t) = 9(t) n e -^P 2 +-o)t $ n (xix; L)$ n (x 2 r, Ije-^™^ 2 '. (82) 

n=l 

Using the results of Appendix [U] we can write 

R (0 \p,x 1± ,x 2± ,t) = e(t)^je-^ 2+ ^^^^ (83) 

with the mean-field expression of the scaling function ^ given in Eq. (|186|) . which, within MFT, 
does not depend on p or £. Within the same approximation the scaling variable (i/To)(£o/ L) z 
is given by 

t= {t/T )(£ /Ly = Qt/L 2 . (84) 

In favor of a compact notation in the following we use this abbreviation keeping in mind that it 
can be replaced by the r.h.s. of Eq. (|84j) . The scaling properties of the response function clearly 
emerge from this expression. Comparing with the general scaling form Eq. (|43|) it is easy to see 
that (p = pL, Xi± = Xi±/L, L = L/£) 

=# (0) (p,xi ± ,x 2± ,t,L) = 2e-(P 2+Z2 )%(x 1± ,x 2± ,t) (85) 

where we used Eq. ()57|) , Moreover, one can easily recover the result for the semi- infinite geometry. 
Indeed, using Eq. (|188j) we find that for Xi± <C 1 (i.e., close to the near wall at x± = 0) and 
i 1 (so that the influence from the wall at x± = L can be neglected near x± = 0) 



(p, xix < L, x 2 ± < L, t) 
9(t) 



To L L 



in agreement with Eqs. (11.18) and (11.19) in Ref. [7]. Note that the previous expression is 
indeed independent of L, as expected for the limit we are considering. Equation (J188|) provides 
also a representation of the response function in terms of the bulk response to a set of image 
excitations. Let us recall that the response function in the bulk is given by 

<J k (p,* 1± ,Z 2± ,t) = ^(t)|l e -(P 2 ^^/ ? 2 ) t /To^^ e -(, 1 x-, 2 x)V(^*/To) (8?) 

To ^/Ant/To 
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x± = x± — L 



x_ 




t f 

confining walls 



X 



+ 



Figure 1: Set of images corresponding to a perturbation actually applied at x\. x 2 is the point 
at which the effects of the perturbation are observed. Image sources are obtained by successive 
reflections of the real and image sources with respect to the confining walls. At each reflection 
the sign of the contribution to the total response changes, starting from a positive sign for the 
actual source. 

According to Eqs. (|83j) and (|188|) one has 

+00 

R {0) (p, x 1± ,x 2± ,t) = Kik(P> <±^2±,t) - <? lk (p, x- ± ,x 2± ,t)} (88) 

n=— 00 

where x^ ± = x\± + 2nL and x~ ± = —x\± + 2nL. The set {x^j_}n^o represents the positions at 
which the "positive" images are located, whereas {s~j_} n gives the positions of the "negative" 
ones. This construction is illustrated in Fig. ^ By virtue of the FDT Eq. I)88|) is also valid for 
the correlation function, with R replaced by C in the expression. Moreover it is an extention 
to the dynamics of the analogous formula known for the static correlation function (see, e.g., 
Subsec. IVA of Ref. [6] and Ref. [37]). 

Let us consider the long-time limit of the response function in Eq. ()82|) . For t > 1, i.e., when 
the effect of confinement is no longer negligible, the sum in Eq. (JH2J) is dominated by the lowest 
mode of the system, i.e., by n = 1. Thus one has 

R {0) (p,x 1± ,x 2± ,t - 00) = e{t)$e-^ 2+L2 y t \^ 1 {x 1 r,m 1 (x 2 r,L)e-^ t + O^ 1 ) 

Jo L 

i.e., the linear response to an external perturbation decays in time exponentially with a factor 
exp[— (q| + L 2 )t\ where, according to our notation, q„ = (p, ttu/L). This exponential decay 
also holds for the two-point correlation functions. Even at T c ^ and for excitations which do not 
break the translational invariance in lateral direction xn, i.e., for p = 0, there is an exponential 
decay due to |qi| > tt/L. This is the result of the combined effect of confinement and Dirichlet 
boundary conditions which suppress homogeneous modes (with vanishing total momentum qo = 
0) in the system. In the case of periodic and Neumann-Neumann boundary conditions for 
this confined system we expect an algebraic decay as function of time for the critical response 
to an external perturbation with p = 0, given that the mode with qo = does occur in 
the corresponding spectra. On the other hand an algebraic decay can be recovered also in 
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the case we are considering. Indeed for T < T c ^ the dimensionless variable L 2 = (L/£) 2 
appearing in the previous equations has to be replaced by — l/2(L/£ t ) 2 within MFT leading to 
the decay ~ exp[— (q 2 — L 2 /2)t\. Accordingly, upon decreasing the temperature T a critical 
value T C {L) < T c ^ exist for which this exponent vanishes. For T = T C (L) the bulk correlation 
length £ attains the value £ c = L/y/2'K. This corresponds to the critical-point shift in film 
geometry, that can also be determined from the onset of a non-trivial order-parameter profile 
(see AppendixEJ). From this point of view the fact that at T = T c ^ the response and correlation 
functions decay exponentially reflects that T c & is above the critical temperature of the system in 
the film, i.e., located within the disordered phase of the film. In the case of a semi- infinite system 
the asymptotic decay of the response function is indeed algebraic even for Dirichlet boundary 
conditions, as one can see directly from Eq. (|86|) for p = 0, T = T c ^, and t/To 3> x 2 ± /^ 2 : 

A(p = 0,x u .,x 2± ,t^oo) = & f±) ^ ^L^L [1 + o(( Xi J^)\T /t))]. (90) 
1 V4vrio \J-qJ so ?o 

This expression is in agreement with the general scaling form in Eq. (|62j) with the mean-field 
values of the exponents and amplitudes (see Eqs. (|38|) and (|57|l). and the universal constant T> 
(see Eq. (|61|l) takes the value 

P<°> = -L . (91) 
v 7r 

This can also be interpreted as reflecting the fact that in the semi-infinite geometry there is no 
critical point shift. The same conclusion can be reached in the case of the film geometry with 
periodic or Neumann-Neumann boundary conditions which, as mentioned previously, within 
mean-field theory do not lead to a critical point shift in films. 

From Eq. (|89|) we can see that for asymptotically large times, the spatial dependence of the 
response function in the film geometry is given by sin(-7rxij_/-L) sin(irx 2 ±/L). 

According to Eq. (|llj) the linear response function R represent the order-parameter profile 
due to a <5-like perturbation applied at an early time. Of course, being derived in linear ap- 
proximation, this function is useful only as long as the subsequent values assumed by the order 
parameter are small enough compared to nonlinear terms. The case of nonlinear relaxation will 
be discussed in Sec. |5] 

In Fig. |21 the function ^(xu_, x 2 _\_, i ), i.e., R^°\p,xu_,x 2 ±,t) up to a spatially constant 
prefactor (see Eqs. (|83|). and (|85|)). is shown for two values of xi±, i.e., the point at which 
the perturbation has been applied at time t = 0. In accordance with our previous observation, 
the response function for t > 1 turns into a sine function with period x 2 ± = 2. We observe 
clearly that there is a qualitative change in the shape of the responding order-parameter profile as 
time increases. In particular the inflection points, which are present just after the perturbation 
has been applied 4 , disappear in the long-time limit, after having reached the closest surfaces. 

In Fig.lHlwe report the time ij(xi±) at which the inflection point (of \P as a function of x 2 ±) 
close to the wall at x± = disappears. Given the symmetry of the problem the analogous time 
for the inflection close to the wall at x± = 1 is simply given by tj(l — xi±). The behavior of 
ti(x~i± — ► 0) can be easily predicted by taking into account that in the semi-infinite geometry, 
for a given x±±, we expect a finite non-zero value tj. According to this argument, in the limit 
L — » oo with fixed the relation 

tl = F tl (x 1A _) (92) 

4 According to Eq. 11881 . for i « 1, ^>(xi±,x 2 ±,i) has a Gaussian form as function of x 2 ± for fixed Six and 
vice versa. 



18 




Figure 2: Time evolution of the scaling function ^(x\±, X2±, t ) which enters into the expression 
of the response function in Eq. (f%5|) (see also Eqs. (|l3*)) and with Xi± = Xij_/L and i = 

(t/To)(£o/ L) z . In (a) the relaxation follows an excitation at the point x\±_ = 0.2 and in (b) at 
the point = 0.5. Reduced times i listed in (a) and (b) refer to the various curves shown 
from top to bottom. 




x 1± 

Figure 3: Time ti(xi±) at which the inflection in the response function disappears as a function 
of the position x~i± at which the perturbation is applied. The dashed curve i(x~i±) = x1 ± /6 
is the quadratic behavior expected for small x~\±_ which actually describes the curve even for a 
wide range of values of x\±. 
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should become independent of L, i.e., (see Eq. (|84j)) 



F tl (y^o)~y 2 (93) 

and thus (see Fig. |3J) 

i T (x 1± -► 0) ~ z 2 ± . (94) 

As discussed in Appendix[E]it is possible to determine analytically, within MFT, the proportion- 
ality factor in Eq. (|94|) . which turns out to be 1/6. Moreover Ft I {y = 1) can also be determined. 

Interestingly, within mean-field theory it is possible to prove (c.f., Subsec. FO)) that ti(x~i±) 
has also the meaning of being the time at which the fluctuation-induced force acting on the 
confining walls is maximal. 

Figures a) and b) clearly show that for t small enough the order parameter profile is 
well localized around the point at which the perturbation has been applied, in accordance with 
the expectation that the effects of the perturbation reach the different points of the system only 
with a certain delay. On the other hand, irrespective of how small t is, the order parameter is 
non-zero in the whole range < X2± < 1 (even everywhere in the case of unbounded geometries), 
as one realizes from Eqs. (|85[l and (j!86[) . This absence of a finite front propagation speed is a 
consequence of the coarse-grained description underlying the field-theoretical approach, in which 
the microscopic time and length scales are assumed to be negligible compared to the mesoscopic 
ones, to the effect that the microscopic dynamics, which of course exhibits a speed limit for 
the front propagation, appears to be actually arbitrarily fast. (This is analogous to the case of 
random- walk models of free diffusion and their corresponding continuum descriptions.) 

For studies of the dynamical properties in the film geometry by means of elastic scattering ex- 
periments one is interested in the two-point correlation function in the (p, x±, (^-representation. 
(The corresponding static properties of thin films near continuous phase transitions have been 
studied theoretically in Ref. [15].) By applying the FDT we can compute the corresponding func- 
tion once the expression for the response function is known in the same representation. In doing 
so we can take advantage of the analytical results known for the static correlation function (see 
Ref. [44] and Appendix[Bj) given that, apart from a factor f2 _1 , the expression for R,( ' in Eq. (|23f) 
is related to that for the static correlation function (i.e., /dw/(27r)C^ (q, uj) = l/(q 2 + <!;~ 2 )) 
by means of a formal shift £ -2 ^ £~ 2 — iuj/Q. Using Eq. H182|) we find 5 

(0) sinh(ax</L) sinh[a(L - x>)] 

aLsmh(aL) (95) 
a 2 = p 2 +£~ 2 - iu/Q. , 

where x^_ = max{xu, X2±} = (x±± + X2±_ + — X2±\)/2, and = min{xij_, X2±} = (xi± + 
X2± — \%i± — ^2_l|)/2. Equation (|93|) agrees with what was found in Ref. [34], Eq. (7), and 
Ref. [38], Eq. (13). Using Eq. (|57|) one can write this expression in the scaling form given in 
Eq. 

/ n \ , sinh(axf ) sinh[a(l — x^)} , , 

^ {0) (p,x 1± ,x 2 ±,0,L =2 1 ±J _ . , L _ l ^, (96) 

a sinn a 

where x ± ' = x ± ' /L, 

a 2 = p 2 + L 2 - id , p = P L, L = L/£, and Q = (uT )(L/t ) z . (97) 



J Note that fixing the branch of the square root defining a = \/p 2 + £~ 2 — i lo/Q, is irrelevant, because Eq. I181H 
is symmetric in a. 
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4.3 Casimir force 



In a confined system the spectrum of the allowed critical fluctuations of the order parameter 
is modified compared to the bulk case, depending on the specific boundary conditions (i.e., 
surface universality classes) and on the film thickness L. This leads to a finite-size contribution 
to the free energy. Accordingly, by varying L one observes that the confining walls are subject 
to an L-dependent effective force F per cross-section area of the film and per It-gT^ (where 
k-Q is the Boltzmann constant) which is the statistical analogue [10, 13] of the Casimir force of 
quantum electrodynamics. In the case of the film geometry with the confining plates belonging to 
(a, 6)-surface universality classes [the case we are currently interested in is the ordinary-ordinary 
(O, O) one] and in the static case F has been shown to scale as function of the thermodynamic 
parameters (up to contributions from additive renormalizations) as [12] 

F(r,h,L) = j ? ^f((h/t ) o)(L/^f^,L/0 (98) 

where the universal exponent 5 and the nonuniversal amplitude f)o are defined through the 
critical equation of state 

m 



h = f) — 
tn 



6-1 

(99) 
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which express the dependence of the bulk order parameter m on the external field h along the 
bulk critical isotherm T = T c ^. Within MFT one has 5 = 3. Because of the two scale-factor 
universality (see, e.g, Ref. [41]) the nonuniversal amplitude rjo is related to the nonuniversal 
amplitudes mo and £o introduced in Sec. [21 via 

fro = Rxi ^ ) (ff )-V , (100) 

tic 

where R c , Kg [see also Eq. (|47|) ]. and R x are universal amplitude ratios whose actual definitions 

and values can be found in Ref. [41]. In Eq. (|98j) the function is a universal scaling 

function with ^ s ^(0,0) = (d — l)A a ^. The amplitude A a ^ is the so-called universal Casimir 
amplitude corresponding to (a, b) surface universality classes. The effective force F between the 
confining walls is attractive when F < and repulsive otherwise. In this sense F can be viewed 
as a special case of the more general case of fluctuation-induced effective interaction between 
different objects immersed in a critical medium. In the following we consider only the case 
(O, O) of ordinary-ordinary boundary conditions and therefore we shall omit the specification 
(a, b) from the scaling functions and amplitudes. Within the Gaussian approximation A is given 
by [11] 

* = -J^ NT (d/2)C(d) (101) 

where T(z) and Q{z) denote the Gamma function and Riemann's zeta function, respectively. The 
universal amplitude A and the universal scaling function j£"( st ) can be determined by computing 
the singular part of the free energy of the system in confined geometry, based on the Hamil- 
tonian JHJ). An alternative approach, more suited for extensions to dynamics, is based on the 
connection between the Casimir force and the expectation value of the stress-tensor T^. (We 
refer the reader to the literature [13] for details.) In particular one finds that the force density 
per kRT c b on one of the confining walls dV is given by 

F=(T ±± )\ 9V (102) 
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where T±± is suitably denned in terms of the order parameter field (p. For the case we are 
interested in one has, in the absence of surface fields (see, e.g., Ref. [13]), 



F i(r)(x\\) = ( T ±±)\dv 




(103) 



Note that in the static case, for which the stress-tensor is a conserved quantity (i.e., d(T^ u ) /dx^ = 
0, where /i, v = _L, ||), {T±±) is actually independent of the point of evaluation, including the 
surfaces. In the dynamic case, in general the force density on the left wall (I) is different from 
the one on the right wall (r) and it may vary spatially along the walls. The mean force per area 
on the wall l(r) is given by 



The relation between the thermodynamic Casimir force (defined from the finite-size behavior of 
the free energy, see, e.g., Ref. [13]) and the expectation value of the stress tensor is based on 
the fact that the equilibrium distribution function of the order parameter field ip is proportional 
to exp{— 7i[ip}}. In the case we are interested in, 7i[(p] is the Landau-Ginzburg Hamiltonian 
in Eq. ([8]). (In turn, the specific expression of the stress-tensor in terms of ip (Eq. (|103j0 is 
determined by TC.) When studying critical dynamics and the effects of time-dependent exter- 
nal fields, such a connection is no longer evident because equal-time correlation functions are 
generated through the distribution oc exp{— W [</?]} only asymptotically for large times, i.e., long 
after any perturbation has been switched off. In principle it is not even obvious how to define 
a thermodynamic Casimir force when studying dynamics, because strictly speaking in this case 
the equilibrium free energy loses its significance. Therefore we assume as the definition of the 
dynamic Casimir force the time-dependent expectation value of the equilibrium stress-tensor, 
which renders the static Casimir force in thermal equilibrium (see also Refs. [37,38]). Heuristi- 
cally, this amounts to assume that at each time, analogous to thermal equilibrium, there is an 
"energy cost" ASL(T±±) associated with the displacement 5L of one of the confining walls, i.e., 
(Tj_j_) provides the local pressure [37,38]. It is desirable to establish a clearer connection between 
this definition of the dynamic Casimir force and the force that can be measured directly in actual 
experiments and molecular dynamics simulations. The previous definition is particularly suited 
for field-theoretical analysis. On the other hand it does not lend itself straightforwardly for the 
study of the dynamic force via Monte Carlo simulations. First, the explicit expression of the 
stress tensor in terms of the order parameter field can be determined in general only pertur- 
batively in terms of the coupling constant go. Second, one has to construct the lattice version 
of the stress tensor in terms of the microscopic degrees of freedom (e.g., spin variables); this 
construction is in general not free from ambiguities. An alternative approach to this problem 
has been taken in Ref. [46] to study the statistical fluctuations of the Casimir force [47] via 
Monte Carlo simulations. However, so far this approach could be implemented only for periodic 
boundary conditions. 

We now consider the case in which the film, thermodynamically close to T c h, is perturbed 
by a time-dependent external field h(x,t). For the ensuing dynamic force density per ksT c £ 
exerted on one of the confining walls one expects a scaling behavior, as in Eq. (|98|) . 




(104) 



Fi(r)(x\\,T,L,t,{h(x,t)}) = 



^^ y) (x||/L,L/e,(i/^^ 



(105) 



22 



in terms of the scaling field h defined via 

h(x,t) = UL/^r Wu h^/L, (t/T )(^/Ly). (106) 

In Eq. (|105j) we have assumed that L does not vary as function of time. For a time- independent, 
spatially homogeneous external field h, reduces to ^^fb introduced in Eq. (|98j). As 

explained in Appendix^! within the Gaussian approximation it is possible to compute the force 
F exactly for a general external applied field h(x,t). In order to elucidate some qualitative 
features of the dynamics of the Casimir force after the perturbation, we consider the particular 
case in which the field is instantaneously applied at a given time t\ and then immediately 
switched off again, i.e., h(x,t) = h(x)S(t — t±) with h(x) localized in the interior of the film. 
After the perturbation the response starts to propagate in the film. At the very early stage the 
response has practically not yet reached the confining walls, so that the force acting on them is 
basically the equilibrium one corresponding to a vanishing magnetic field. In course of time the 
perturbation induced by the field hits the confining walls and correspondingly the force exerted 
on them changes. Finally, because of the relaxational character of the dynamics, in the limit of 
long times the perturbation induced by h disappears and accordingly the effective force reaches 
again its equilibrium value. 



4.3.1 Planar perturbation 

In order to illustrate such a behavior we consider the case in which the perturbation does 
not break the translational symmetry along the confining walls and is localized in the plane 
x± = xi±, i.e., h(x,t) = h\y5(x± — xi±)S(t — t±). From Eqs. (|231|) and (|232|) one finds for the 
left wall (the upper index (°) indicates the Gaussian approximation) 

^ dym (L,i,h w ) = ^°\o,L) + h 2 w \d^ ± ^ \p = o,x 1± ,x 2± ,t-t 1 ,L)\ S2±=0 \ 2 , (107) 

o L 

where the scaling function of the response function is given by Eq. ()85|) and 

Pi Arh w (108) 

so/ 4Vo 

is the scaling variable associated with hyy. In the present case the functional dependence of the 
force density on h(x,t) reduces to the dependences on hw, xi±, and t\ and there is no spatial 
dependence of the force on the lateral coordinates. According to Eq. ()1U7[) . the response of the 
Casimir force to an external field is related to the square of the spatial derivative of the response 
function evaluated at one of the confining walls. As expected, the Casimir force depends only 
on the time St = i — i\ elapsed since the application of the external perturbation. 

In the following we discuss in more detail the relaxation of the Casimir force at the bulk criti- 
cal point T = T Cj b, i.e., for L = 0. In this case &(°\p = 0, x\±, x 2 ±, Si, L = 0) = 2^f(xxx, x 2 j_, Si) 
[see Eq. ()85(l ] and plots of the function ^ are provided in Fig.^lfor various values of its scaling ar- 
guments. In the present context the square of the derivative of these graphs at x 2 ± = matters; 
x~i± is the position of the applied perturbation. Moreover, the static contribution in Eq. (|107|) 
is simply given by the Casimir amplitude A [see Eq. IjlOlJ) ]. i.e., ^ st )(°)(0, L = 0) = (d — 1)A. 
Note that A < for the case we are interested in [11], whereas the contribution stemming from 
the field perturbation is always positive. Therefore the overall sign of the effective force on the 
left wall depends on the strength hyy of the applied field and on the actual position where it has 
been applied; in particular it may become positive, i.e., repulsive within a time window during 
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the relaxation process. Given the expression of \& in Appendix [(""] [see Eq. (|186|) ] it is possible 
to compute analytically the critical scaling function in Eq. ()107j) . According to the qualitative 
behavior described above and according to Fig. [""J \d X2± ^°^(p = 0,xi±,X2±,i,L = 0)\ X2±=0 \ = 
2\dx 2± ^(x~i±, x~2_l, i )\x 2± =o\ displays, as function of time, a maximum for i = tM(x~i±) which is 
implicitly defined by the condition 

dtdx 2± -$(x 1± , x 2 ±, i )U 2X =o,t=f M (Six) = • ( 109 ) 
On the other hand, from the definition of ^ in Eq. 1)183)) . together with Eq. ()18|) it is easy to real- 
ize that dt^(x~i±,X2_L,i) =_di 1± ^(xi_L,x 2 ±,i) = d^ 2± ^(xi ±,x 2 ±,i). Accordingly dtd X2± ^(x±±, 
%2±,t) = 5| 2± ^ / (xij_, x~2±, t ) and the necessary condition ()109)) can be written as 

dl 2x y{x lx ,X2x,t)\x 2±=0 ,i=t M ($ 1± ) = ^ 3) (xi±,0,t M (x 1± )) = 0. (110) 

where we use the notation introduced after Eq. 1)215)1 in Appendix[El Comparing Eq. (jllOj) with 
Eq. (|221|) one sees that tM^ijJ and (Eq. (19*2*)) ) satisfy the same equation. On the basis of 

the qualitative behavior of ^ one expects the solution to be unique and thus ^m(^i±) = ti(x\±), 
i.e., the time at which \d X2± ^(xi±, x~2x, i )\ X2±=0 \ displays a maximum for a fixed value of X\± 
equals the time ii(x~i±) at which the inflection point of ^(xi±, x~2±, t ) as a function of X2± 
reaches the surface at X2± = (see Fig. |3J) . 
Motivated by Eq. ()1U7|) we define 

lWf . ,_^ (dy)(0) (o,i/(^),W)-^ (st)(0) (o,o) 



hi 



1 _ o 

- [d X2± ^(x 1± ,X2±,t I (xi±))\ X2±= o] 



w (111) 



^(xi±,t) = ^ — ToTJ ~V — — ( 112 ) 



so that (d—l)A + hyyA^(xi±) is the maximum value of the force to which the left wall at x± = 
is subject if the field is applied at the normal distance x\±. Figure [3] shows the dependence of 
the maximum of the field-induced force on the position x\± of the perturbation whereas in Fig. [3 
we show the time dependence of the normalized dynamic part of the force 

^ (dy)(0) ( o,i,h w )-& ww (o, o) 

h 2 w A^(x 1± ) 

for various values of x\±. The asymptotic behaviors of A^(x\±) for x±± — > and x\± — ► 1 
are given by A%(x 1± -»■ 0) = Ki/xf ± with K\ = 27e~ 3 / 7r ~ 0.427889 and A%(x 1± -»■ 1) = 
ax(l — xi±) 2 with a\ ~ 17.541, respectively (see Appendix IF.2() . As Fig. [3] clearly shows, 
they provide very good approximations to the actual function already for x\± < 0.6 and 
xi± > 0.75, respectively. As expected, when the plane in which the external field is applied 
approaches the distant wall (i.e., x\± — ► 1), the actual response of the system is reduced due to 
the Dirichlet boundary condition there and, due to the relaxational character of the dynamics, 
the ensuing perturbation affects only slightly the wall at x± = 0. This qualitative behavior 
clearly emerges from Fig. |3J On the other hand, if the external field is applied close to that wall 
where the response is monitored, the spatial variation of the response function as a function of 
X2_l is sufficiently pronounced to induce a strong force on the close wall. However, this force 
quickly decays with time, and indeed ii(xi± — > 0) — » (see Fig. H"}. The time dependence of 
the normalized dynamic part of the force (x\_\_,i ), reported in Fig. [""J depends only weakly 
on x\± when plotted as a function of the scaled time i/ii{x\±). In Eqs. 1)245)1 and "246)1 we 
provide the asymptotic expressions for ^(i u — > 1, i) and ^(ijx — ► 0,t), respectively, 
shown in Fig. ["0 For x\± fixed and within mean-field theory decays asymptotically as 
^{xi±,i^ oo) ~ e~ 2 ^ 1 (see Eq. $£Dj ). 
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Figure 4: Dependence of the amplitude AW(xi±), which determines the critical Casimir force 
maximum (d— l)A + hy V A^(xi±) on the left wall [Eq. QUI)) ], on the normal distance x\± where 
the planar perturbation is applied. The asymptotic behaviors for x±± — > and x±± — > 1 are 
indicated as dashed and dashed-dotted lines, with K,\ = 27e~ 3 /7r ~ 0.427889 and a\ ~ 17.541, 
respectively. Note that already for x\±_ < 0.6 the asymptotic behavior IC\/xf ± for xi± — ► 
provides a rather good estimate of the actual value A^(xi±). 




8t/ti(x\±) 



Figure 5: Time dependence of the normalized dynamic part of the critical Casimir force 
&g(xi±,i) = [^ dym (0,i,h w ) - ^( st )(°)(0,0)]/[/i^^(xi ± )] for various values of the po- 
sition xi±_ where the perturbation is applied. For t = ii(xi±) the Casimir force reaches its 
maximum (compare Fig. EJ). The dashed and dash-dotted lines indicate the limiting shapes of 
the curves for x\± — > and x\±_ — > 1, respectively. For i S> l/(37r 2 ) the actual curves decay as 

c -2n 2 t 
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4.3.2 Localized perturbation 

Here we discuss the case in which the external field is localized at a single point (xiii , x\±) within 
the film, i.e., h(x, t) = hpd(x\\ — x^)5(x± — xi±)5(t — ti). From Eqs. (|231j) and (J232() one finds 
for the left wall at x± = 



jr( S t)(0) (0;Z) + [d S2 ^(x 1± ,x 2± ,Si)\ S2±=0 ] 2 



-{S5i ll ) 2 /(2St)-2L 2 St (113) 



where St = t — t\ and Sxi\ = xn — x^. The function is given by Eq. (|186|) and 

#1 "^-^ • (H4) 



hp 



In the following we discuss in more detail the relaxation of the Casimir force at bulk criticality 
T = T C] b, i.e., L = 0. Different from the case of a planar perturbation, here at a given time 
the force varies laterally. It depends on the normal distance of the epicenter from the wall 
under consideration and on the radial distance |<5x||| from the epicenter projected on this wall. 
The force decreases monotonically for increasing radial distances \5xn\. The qualitative time 
dependence of the force, generated by a point-like perturbation, is expected to be independent 
of the actual lateral position where it acts: As in the case of planar perturbations, the force 
equals the equilibrium one for very short and very long times with a maximum in between at 
t = Im (£xm , x\±) measured from t\. Upon increasing the lateral distance \8x\\\ from the source 
of the perturbation Im is expected to increase, because the perturbation has to cover a larger 
distance until it hits the wall at the specified point. The asymptotic behavior of Im f° r large 
\Sxu \ can be inferred from Eq. ()113|) by taking into account that 

oo 

^ (0,1 W,0,t") = d S2± y(x l± ,x 2± ,i)\ S2±=0 = 27r^ne-^™ 2 *sin(vrnS 1± ) . (115) 

n=l 

For 7r 2 t > 1 the sum is dominated by its first term. This allows one to determine its extremum 
as function of t. At leading order the result is independent of xu_: 



4vr 2 (c5x||) 2 + (cZ-l) 2 -(ci-l) 



tA/(|<5x|j | > 1, xi±) = - 1 — 2 . (116) 

It is remarkable that in spite of the diffusive propagation of the perturbation (Eq. (|113|0 . the 
maximum of the laterally varying force moves asymptotically with constant velocity: tj\/(|(5x|| | — > 
oo) = \Sxu | / (2tt) so that with Eq. (|84f) this asymptotic speed v is given by 

v = 2<n{^/T ){^/L) z - 1 . (117) 

Thus this speed decreases with increasing film thickness. It is possible to provide an estimate 
of the typical value of v by considering the values of £o an d To that have been experimentally 
determined. In Ref. [48] the critical dynamics of a ultrathin film (bilayer) of iron grown on a 
tungsten substrate has been investigated via the magnetic ac susceptibility. This system should 
provide a realization of the two-dimensional Ising universality class with Model A dynamics 
(z ~ 2.1 [40,48]). In particular To~ exp for the exponential relaxation time (see Appendix lA"|) 
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has been obtained from the fit of experimental data, yielding Tq =2.6±0.6x10 s. For 
£o we take £o ~ 3A corresponding to the lattice constant, as it is usually the case in magnetic 
materials [39]. Accordingly, for a thin film with L = 50 £o> one finds v ~ 0.1 m/s. 

By using the behavior of ^I^ ' 1 ) for ir 2 i <C 1 reported in Eqs. ()236[) and (|237[) for x±± — > 1 
and xi± — > 0, respectively, one can determine the corresponding behaviors of tyi under the 
assumption n 2 iM *C 1 for which they are given by 



d + 5 + (ix,,) 2 - . [d+ 5 + (Sin) 2 ] 2 - 4(d + 2)[(ax, )2 + 1] 

f » = " WT2) (118) 

for x~i± — > 1 and 

((5X||) 2 + X 2 

t M = 1 — !ii -ii f or ^ _> o . 119 

2(d + 2) v 7 

Comparing with the numerical determination of tjv/(^ x | , %i±) it turns out that Eq. (|119|) provides 
a good approximation for tM up to few percent in the region |<5x|||,xij_ < 0.5, which corresponds 
to £m < 0.1. (Note that for d = 1 and Jxii = we formally recover the result for the planar 
perturbation, as can be seen by comparing Eqs. (|1()7|) and (|113[) .) In particular, from Eq. (|118j) 
one finds i M {^\\ = 0,lu 1) = (4 - V / TT)/10 ~ 0.0683375 for d = 3 and tM^xy = 0,xi ± -> 
1) = (9 — \/57)/20 ~ 0.0604236 for c? = 4 for which the Gaussian approximation becomes exact 
apart from logarithmic corrections. In the case x\± — > one finds from Eq. (|119|) . ^(Jxii = 
0,x\± — ► 0) = x 2 _|_/[2(<i + 2)]. Figure a) shows ijw (5xn , in) in d = 3 as a function of \65L\\\ 
for certain values of x±±. The asymptotic behaviors for x\±_ — > 1 (Eq. (|118jl ) and for x\±_ — > 
(Eq. (|119|l ) as functions of |5xm| are also shown as dashed lines up to a corresponding value 
of t,M — 0.1 (fulfilling the condition 7r 2 iM ^ 1 under which Eqs. (|118j) and (|119|) have been 
derived). The dash-dotted line in Fig. Efa) indicates the leading linear behavior of tM for large 
|<5x||| [see Eq. (|116j) ]. which provides a rather good approximation of the actual dependence 
already for |5xm| > 1.0. Figure|H{b) shows the comparison between Fm^Xii = 0, xi±) in d = 3,4 
and ij(xi±) (see also Fig. |3J) [formally corresponding to the case d = 1]. The dashed curves 
indicate the approximate expression Eq. (|119|) valid for small x%±, actually providing a rather 
good approximation already for x\±_ < 0.6 — 0.8. Figure HJb) clearly indicates that tM decreases 
upon increasing d. Indeed, for fixed x±±, the factor in the denominator of Eq. (|113j) increases 
the force at short times compared to the case of the planar perturbation (Eq. (|107j0 where it is 
absent. This causes the maximum of the resulting force to occur at earlier times i, as indicated 
by Fig.Etb). 

In analogy to Eq. (|llll) . we define as a measure of the maximum force 



^ (dy)(0) (0, t M (£x|| , x 1± ),h P ) - (0, 0) 



A%(5x\\,xx±) 

>bp (120) 

l e -(«*||) 2 /[2fAf(**||,Sl±)] 2 



2 [47rt M (5x||,xi ± )] f 



— [dx a j_^(xi±, x 2 ±, t M (S*\\ , Xl±))\x 2± =o] 



We first discuss A^(5xii = 0, in), i.e., the maximum of the field-induced force at the point of 
the wall closest to the source of the perturbation. Figure [7| shows the function Aa(<5x|| = 0, x\±) 
for d = 3, together with its asymptotic behaviors determined in Appendix 11' .21 According to 



Eqs. (ESOJ) and (t23T|) they are given by A%(8x\\ = 0,x 1± -► 0) = /C d /x 2 ^ +1) , with K d defined 
in Eq. (J2HJ), and = 0,x ix 1) = o d (l - xi_l) 2 , with a d defined by Eq. (1232]) . The 

qualitative dependence of A£(5xm = 0,x\±) on x\± is the same as that of A^(xi±). 
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Xxx. 

Figure 6: Time iM(<^||,%j_), at which the critical Casimir force at the point on the left wall 
at X]M + <5x|| reaches its maximum, as a function of |Jxn| in d = 3 (a). The external field has 
been applied at (xx\\,xx±). The solid lines, from top to bottom, refer to the cases x±± = 3/4, 
1/2, and 1/4, whereas the dashed lines show the asymptotic behaviors for x±± — ► (Eq. (|119|l ) 
and — > 1 (Eq. I|118|) ) for 7r 2 tM <C 1. The dash-dotted line indicates the approximation of Im for 
large |<5xn| given by Eq. (|116|) . which is already rather good for |5x||| > 1.0. For large |5xm| the 
time ?m becomes independent of x\±] it increases linearly with the slope giving the inverse of 
the speed for the propagation of the force maximum on the left wall (see Eq. (|117|l ). In (b) we 
compare the time Im{^\\ = 0, xuj for d = 3, 4 with tj{x\x_) (see Fig.[3J), formally corresponding 
to the case d = 1. The dashed curves indicate the quadratic behaviors expected for small xxx. 
and arbitrary d (see Eq. (|119(l ) which actually describe the curves rather well even for a wide 
range of values of x~xx_- 
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Figure 7: Dependence of the amplitude A^(#Xn = 0, x\±), which determines the critical Casimir 
force maximum (d — 1)A + hpA^(65L^ = 0, at the point of the left wall closest to the point 
within the film where the external field was applied, on the normal distance of the perturbation 
from the left wall. The curves refer to the three-dimensional case, although the qualitative behav- 
ior is the same in d = 4. The asymptotic behaviors for x\_\_ — > and x±± — > 1 (see Appendix lF.2|) 
are indicated as dashed and dashed-dotted lines, respectively, with /C3 = 5 5 /(47r 2 e 5 ) ~ 0.169773 
and 03 ~ 17.927. Note that, as in the case in Fig. 01 the asymptotic behavior K,^/x\^ for 
xi± — > provides a rather good estimate of the actual values already for x%± < 0.65. 

According to Eq. (|113j) . at a given time 8t elapsed after the perturbation, the field-induced 
part of the Casimir force is a decreasing function of |<5xm|. Therefore the maximum amplitude 
vl^(<5x|| , xi^) decreases with increasing lateral distance |<5x||| from the point where the field has 
been applied. Figure |S] shows the dependence of the maximum force amplitude ^4£(5x|| , x\±) 
in d = 3 on <5x|| for fixed values of normalized to its value at 5xii = (compare Fig. EJ). 
The solid lines correspond, from top to bottom, to x\±_ = 1, 3/4, 1/2, and 1/4. For xu_, |<5x||| < 
0.5 the approximate expression of j4£(<5x||,xi^)/A£(c)X|| = 0,xu_) ~ [1 + (<5xii ) 2 /x 2 _|_]~( d+2 ) 
(Eq. Q254j0 provides a very good approximation of the actual curves. For |<5x||| — > 00 this ratio 
decays ~ exp(— 27r|5x|| |) (Eq. (12531) ). In Fig.|H]the dashed line refers to this asymptotic behavior 
for x\±_ — > 1. It is interesting to note that for all x\± the maximum of the field-induced force 
decays rapidly with \S5tu | and is practically negligible compared with its value at <5x|| = already 
for |(5x|| I ~ 1. 

Figure IH1 shows the time dependence of the normalized dynamic part of the Casimir force 

„ P , ^, (dy)(0) (P,i,h P ) - ^ st )(°)(o,o) , , 

&Z(65t h x 1± ,t) = -J y - // V > (121) 

for d = 3, fixed xu_ = 1/2, and various values of |5x||| = 0, 1, 2, and 4, from top to bottom. 
When the jF^" is plotted as a function of t/tu its shape does neither depend sensitively on the 
specific value of x\± nor on the value of |5xy |. For fixed <5x|| and x±±, decays asymptotically 
as &£(6xpx 1± ,i 00) ~ (5i)-( d -Ve- 2n2sI (see Eqs. (tTT^I and (tTT5l ). 
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Figure 8: Dependence of the normalized amplitude ^4£(<5x|| , xx±) of the critical Casimir force 
maximum with d = 3 on |<Jxn | for various fixed xi±. The dashed line is the asymptotic behavior 
for |<5x||| 3> 1 and x±± — » 1, as obtained from Eq. H253|) , For \6xu\ — > oo the curves vanish as 
~exp(-27r|5x|||) (Eq. (ESSJ)). 
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Figure 9: Time dependence of the normalized dynamic part of the critical Casimir force 
&K(x 1± ,i) = [^ dym (0,i,h P ) - ^( st )(°) (0,0)]/ '[h 2 P A%{55z h x 1± )} for various lateral distances 
\Sx\\ \ from the location of the perturbation at x%± = 1/2, with d = 3. For i= tM(S5c^,xi±) the 
Casimir force reaches its maximum (compare Fig.EJ). The different curves refer to \Sx»\ = 0, 1, 2 
and 4, from top to bottom. The qualitative behavior of the curves is the same for other values 
of xi±. For 5t — > oo the curve vanish as (5t)~^ d ~^ 'e~ 2w 51 . 
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4.4 Correlation function 



4.4.1 General expressions 

As discussed in Subsec. 13,21 the FDT relates the correlation function and the response function 
via Eq. (|31|) . Thus, taking into account Eq. ()83[) . one finds that the correlation function can be 
cast into the scaling form given in Eq. 1)4 5 jl with 

tf {0) (P,x 1± ,x 2 ±,i,L) = 2 dse- s ^ L ^(x 1± ,x 2L ,s) (122) 

J\t\ 

where the value of the non-universal amplitude be is given by Eq. 1)56)1 and where we have 
introduced 

to = pr^p . (123) 

The dependence of the scaling function on t for selected values of x\± is shown in Fig. HOI at 
criticality and for p = 0, i.e., to = oo, and for to < oo in Fig. 1111 As expected, vanishes in 
the limit t — » oo. According to Eqs. 1)182)1 and ()15j) the scaling function ^^\p, x-\j_,X2±,L) = 
<?f (°)(p, x 2 ±, i = 0, L) of the static correlation function is given by 

«*0), - _ o sinh(tV 1/2 x<)sinh[ty 1/2 (l - s>)] 

®st {P,Xi±,x 2 ±,L) - 2 — — ^ 1/2 , (124) 

*o sinh(t ) 

which, independently of the value of to, is characterized by a cusp-like singularity for x\± = x 2 j_, 
clearly displayed in Figs. HUI and HTl In particular at bulk criticality, i.e., L = and for p = 
(i.e., to = oo), one finds (compare Eq. (|95j)) 

^ s ( t 0) (p = 0,x 1± ,x 2± ,L = 0) = 2x<(l-x>) , (125) 

which corresponds to a triangularly shaped correlation function vanishing at the boundaries [see 
Fig. and Eq. 1)182)1 ]. The two-point correlation function in thin films can be probed by scat- 
tering X-rays or neutrons under grazing incidence [15,49]. The dependence of the corresponding 
scattering cross section on the lateral momentum transfer is dominated by the singular behavior 
of the two-point correlation function in planes parallel to the surfaces of the film. With future 
technologies it might be possible to extend such kind of scattering experiment into the time 
domain. Therefore we discuss the particular case x^± = x 2 ±, i.e., C(p, x±, x±, u). The analysis 
of this particular case may also serve to enhance the general physical insight into correlations 
in films. 

In order to understand the behavior of this quantity we shall consider two relevant limits: 
(a) the behavior close to the confining walls, i.e., x±_ <C L (surface behavior; see Appendix [HJ, 
and (b) the behavior in the interior, i.e., x± = L/2. Crossover effects are expected to occur in 
between. 



4.4.2 Behavior near the wall 

To this end we rewrite Eq. 1)95)1 . use the FDT (see Eq. 1)34))) . cast it into the scaling form given 
in Eq. 1)44)) . and use Eq. ()56)) for the nonuniversal amplitude: 

^io) ( - - - — j \ 2 T cosh[a(l- |xi_L-X2±|)]-cosh[a(l-xi_L-x 2 ±)] / 10 «n 

^ u, (p,x 1± ,x 2 ±,uJ,L) = -Im _ . r _ (126) 

u) a sinn a 
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x 1± 1 x 2 ± xi± 1 x 2 ± 

(a) (b) 

Figure 10: Time evolution of the mean-field scaling function (p, x±±, X2±, i, £) which enters 
into the expression of the correlation function in Eq. (jUjl with x~i± = Xi±/L, i = (t/To)(^o/ L) z , 
L = L/£, and to = l/(p 2 + L 2 ) for p = and T = T c ^, i.e., to = oo. We show correlations 
between points X2± and x±± = 0.2 (a) and x±± = 0.5 (b). Reduced times i listed in (a) and (b) 
refer to the various curves shown from top to bottom. At i = 0, the static correlation function 
fft°) exhibits a cusp at = X2±- 
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where a has been defined in Eq. ()97[) . By using the previous equation we obtain 

^roV- - - - f\ 2 x cosho - cosh[o(l - 2x_l)] 

>{p,x±,x±,u),L) = — Im = -— 

uj a smh a 

so that for |a|rEj_ = |a|xj_ <C 1 (i.e., sufficiently close to one of the walls) 

&(°Hp,x±,x±,u;) = -4^Im{acotho(l + 0(\a\x±))} , 

UJ 

in agreement with the general behavior of ^ close to one wall, given in Eq. (|63[) with 



AO) 
7 w 



p, uj, L) = — — Im {a coth a(l + 0(\a\x±))} 



uJ 



(127) 



(128) 



(129) 



We can distinguish two regimes, i.e., \a\ = \a\L 3> 1 (relevant for discussing the connection with 
the results for the semi-infinite geometry) and \a\ = \a\L <C 1. 

• |a| = |a|L^>l:In this case a coth a = a(l + 0(e~ 2a )) (we choose the branch with positive 
real part in the square root defining a, see footnote EJ) and thus 

(PjZL^-U^ 1 ) = -^Im{a(l + 0(|. 

UJ 

Using the definition of a, one obtains 

,(0)/_ - - - r= 1 



(130) 



p 2 + L 2 + ^(p2 + L2)2 +tD 2 



(131) 



■ 4 /(p 2 +L 2 ) 2 +cu 2 » 1 and x± {/ (p 2 + L 2 ) 2 + Q 2 < 1 . 



Equation (|131|) renders the limiting behaviors listed in Eqs. (|69|) . (|70j). and l|71|l with the mean- 
field universal amplitudes 



A W(0) _ 2 
C W(0) _ 2 



(132) 
(133) 
(134) 



and the mean-field values of the critical exponents. Corrections to the formulae Q69J) . (|7U|1. and 
(|71|) are exponentially small. 

• |a| = |a|L <C 1: In this case one has (see Appendix O Eq. (|261|) ). 



^ (0) (P, 



x±,x ± ,u>) =-x\{l - x L ) 2 { 1 - — (1 + 2x ± - 2x 2 ± )(p 2 + L 



15 



(P Z + L 2 ) 2 - -uj* 



+ ... 



{/ {p 2 + L 2 ) 2 + uj 2 « 1 , x ± <l. 
Near one wall, i.e., for xi <C 1 this expression leads to 



(P,X±,X_L,U) 



i _ A ( p2 + m + J_ ( p2 + -j}f _ JLrf + . . . 

15 VP ; 105 VH 7 315 



(135) 



(136) 
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2 4 6 8 10 

y = L/C [y = \p\l] 



Figure 12: Mean-field scaling function ^/(0,0,y) ["^/(y, 0, 0)] which enters 
into the expression of the correlation function C(p = 0,x±,x±,uj = 0) = 
Oc(L/^) l -^ z (x ± /L) 2 ^y^ w (0,0,y L/0 [C txit (p,z ± ,z_ L ,u - 0) = 

o C {L/^) l -^ z {x L /L) 2 ^-^l^ w {y = P L,0,0)] for L/x± » [L/x± » 1, \p\L], so 

that V$(p,0 t y -> 0) [^(y - 0,0,0)] = (4/3) [1 - (2/15)y 2 + (2/105)y 4 + 0(y e )} and 
^ (0,0,1/ -> oo) [^(y -> oo,0,0)] = 2/y + 0(1/?/ 2 ). Beyond mean-field theory the 
asymptotic behavior of "^vk(0, 0, y = L/£) and ^wiy = pL, 0, 0) for large y is given by Eqs. Ijfifijl 
and 1)68(1 , respectively. The index W of the scaling function here and in Fig. ^3 indicates the 
behavior near the wall. 

From this expression we can identify the mean-field value of the universal constant introduced 
in Eq. ((Ml): 

A W(o) = 4/3 _ (137) 

Figures IT2l and fTTTl display the functions (p = 0,lu = 0, Z), "j?^/(p,cD = 0, L = 0), and 
(P = 0,u>, L = 0), respectively, together with their asymptotic behaviors given in Eqs. (|69|) - 
G3J and ftT^jl . 

Thus, in contrast to the case of a semi-infinite geometry (in which the validity of the behaviors 
in Eqs. (|70|l. and (|71j) extends down to l/£ = 0, u; = 0, or p = 0, respectively), the critical 
correlation function in the film does not diverge upon approaching the origin of the (£~ \p,u)- 
space. This is a consequence of the critical point shift in the film geometry. Of course this 
divergence is correctly recovered in the limit L — > oo (i.e., in the semi-infinite limit), where this 
shift is absent. 

4.4.3 Film behavior 

Now we consider the behavior of C^(p, x±, x±, lj) in the middle of the film, i.e., C^°'(p, L/2, L/2, ui). 
According to Eqs. (|44|) and (|77[) this amounts to studying (p, u), L). In that case we find 
(see Eq. RTTH ) 

^ ( °W,Z) =^°)(p, 1/2, 1/2, 0, L) = -Im tanh( " /2) . (138) 

u a 

Again, we consider the two possible regimes \a\ = \a\L ^> 1 and \a\ = \a\L <C 1. 
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y = ujT (L/Co) z 



Figure 13: Mean-field scaling function (0, y, 0) which enters into the expression of the correla- 
tion function C crit (p = 0,x±,x±,u) = d c (L/^) 1 - v ~ z (x±/L) 2 ^-^^ w (0, y = ujT (L/^Y,0) 
for T = T C)b and L/x x » 1, uT (L/^o) z , so that tf${0,y -» 0,0) = (4/3) [1 - (2/315)y 2 +0(y 4 )] 
and«$?(0,y-» oo,0) = 2 v / 2/y 1/2 + 0(l/y). Beyond mean- field theory the asymptotic behavior 
of ^V(0,y = wT (L/^ ) 2 ,0) for large y is given by Eq. (|67|) . 



|a|L 3> 1: In this regime one has 



^°\p,u),L) = ^Im<j 3(1 + 0( 



-25 



)) 



so that 



(139) 



f ) (p,cD,Z) = V2- 



p2 + £2)2 + ^2 > j _ 



(140) 



From this equation one easily recovers the limiting behaviors given in Eqs. (|74|). (|75|). and (|76|). 
with the mean-field universal amplitudes 



= 1 
= ^2 
f - 1 



(141) 
(142) 
(143) 



and with the corresponding mean-field values of the critical exponents. 

• \a\ = \a\L <C 1: On the other hand, in the limit \a\ = \a\L <C 1 one can make use of Eq. (|135|) 
(valid in the case \a\L, \a\x± <C 1 and arbitrary x± < L), with x± = 1/2, leading to 



f>(p,Q,L) 



12 



1-I(p 2 + L 2 ) + — (p 2 + Z 2 ) 2 
5 V ; 560 v ; 



!7 2 



1680 



(p 2 +Z 2 ) 2 +cD 2 « 1 . 



(144) 
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From Eq. (|144|) one infers the universal amplitude 

A 1 ® = i (145) 

introduced in Eq. 1)73)1 . As it is the case near a wall (see the remark at the end of Subsec. 14.4.2)1 . 
also in the middle of the film the correlation function does not diverge at the origin of the 
(p, w)-space; divergences appear only in the limit L — > oo, for which the critical-point shift 
disappears and T c & coincides with the critical temperature of the system. In particular, from 
Eq. (|135|) one finds 

^ c ( 2(p = 0, xl, x Xl Q = 0) = ^1(1 - x ± ) 2 . (146) 

In Fig. ^] we show the functions ffj (p = 0, Q = 0, L) and 'tfj (p, u) = 0, L = 0) which have 
(within MFT) the same scaling function if the scaling variables are appropriately chosen. (In 

Fig. Elthe quantities in square brackets refer to (p,uj = 0,Z = 0).) The limiting behaviors 
of this function as given by Eqs. (|7H) [ (fTo]) ] and (|144|) are indicated as dashed and dash-dotted 
lines. In Fig. 1151 we plot ^j°\p = 0,Q,L = 0) and its limiting behaviors given by Eqs. (|144j) 
and (|75|) . The comparison between the scaling functions near the wall fFigs. WH and fT3* )l and in 
the interior of the film (Figs. [T4l and 115 )1 shows that they exhibit the same qualitative behaviors. 
However, the scaling functions in the interior decay more rapidly for large scaling variables 
than their counterparts near the wall. Moreover, their absolute values scale very differently as 
functions of L. 




V = L/C [y = \p\L] 



Figure 14: Mean-field scaling function (0, 0, y) (y, 0, 0)] which enters into the expression 
of the correlation function C(p = 0,x ± = L/2,x ± = L/2,u = 0) = 6 c (VCo) 1 ^^/(0, 0, y = 
L/0 [C crit (p,x ± = L/2,x ± = L/2,u = 0) = 6 c -(L/e ) 1 ^^(y = p£,0,0)], so that 
^ (0) (0,0,y -> 0) [^ (0) (y - 0,0,0)] = (l/12)[l-(l/5)y 2 + (17/560)y 4 +O(y 6 )] and ^ (0) (0, 0, y -> 

oo) [^j (y — > oo,0,0)] = 1/y 3 + 0(l/y 4 ). Beyond mean-field theory the asymptotic behavior 
of ^/(0, 0, y = L/C) and ^/(y = pL, 0, 0) for large y is given by Eqs. (|74j) and (|76j). respectively. 
The index / of the scaling functions here and in Fig. El indicates the behavior within the interior 
(i.e., middle) of the film. 
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Figure 15: Mean-field scaling function (0, y, 0) which enters into the expression of the corre- 
lation function C crit (p = 0, x± = L/2,x± = L/2,u) = 6c(V^o) 1_J? "^i(0, y = uT (L/£ Q ) z , 0) 
for T = T c>b , so that ^ (0) (0,y -> 0,0) = (1/12) [1 - (17/1680)y 2 + 0(y 4 )] and ^j \o,y -> 
oo,0) = \/2/y 3 / 2 + 0(l/y 2 ). Beyond mean-field theory the asymptotic behavior of ^i(0,y = 
ujTq(L/^q) z , 0) for large y is given by Eq. (|75|). 

5 Nonlinear Behavior 

In the previous sections we investigated the linear response and correlation functions within the 
Gaussian approximation for the dynamical functional (Eq. (|1U|) ). The response function captures 
the behavior of small perturbations around the equilibrium state of the system. The response 
function analyzed in Subsec. [4.21 is useful to describe the relaxation process from an initial state 
with a small value of the order parameter such that nonlinear terms can be neglected. On 
the other hand both in experiments and numerical simulations this is often not the case given 
that it is much easier to investigate the response of the system to finite changes of the control 
parameters (such as temperature and external fields). While typically the former relaxation 
process is characterized by exponential decays in time, the latter leads to power-law behaviors. 
As a concrete example one can think of an experimental protocol in which the system is initially 
prepared in an equilibrium state with a non-zero value of the order parameter. This can be 
realized, for example in magnetic systems, by preparing the sample in its low-temperature 
phase or by applying an external field. Then the external parameters are changed in a way that 
the order parameter in the corresponding equilibrium state vanishes, as it is the case above or 
at the critical temperature, in the absence of external fields. In order to describe the ensuing 
relaxation it is crucial to account for the effects of nonlinear terms, as will be discussed below. 

Let us now consider in more detail the case in which an external field /i(x, t), coupling linearly 
to the order parameter in the Hamiltonian TC, is present during the relaxation. The evolution 
equation for the quantity m(x,t) = (ip(x,t)) can be derived in a standard way [50,51] and, at 
the lowest level in the loop expansion (tree- level), it reads 

h(x, t) = d t m(x, t) + Q[-A + + ^rn 2 (x, t)]m(x, t) . (147) 

In the case of confined geometries this equation, which is valid in the bulk, has to be supple- 
mented with the proper boundary conditions. For bulk systems with homogeneous external fields 
the linear regime is identified as that one in which the nonlinear term in Eq. (|147|) is negligible 
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compared to £ 2 , i.e., ^-m 2 (x, t) <C £ 2 . In view of Eq. (|4Tj) this means m(x, t)/mo -C £o/£- 
(At the ordinary transition, the equilibrium magnetization close to the surface vanishes as 
~ (T C)6 - T)^ 1 for T -> T~ b , with /3 X = 1 within MFT. Therefore, sufficiently close to T C)h , 
the previous inequality is always fulfilled. Nevertheless the spatial inhomogeneity caused by the 
boundary condition at the surface yields a contribution - due to the term —Am in Eq. (|147[) - 
which effectively reduces the one due to the linear term previously considered. Eventually, the 
nonlinear term dominates and the proper nonlinear relaxation is displayed.) In that regime 
and in the absence of an external field the temporal relaxation from an initial state with a 
homogeneous order parameter is exponential. The range of validity of the linear approximation 
shrinks upon approaching the critical point [50-53]. At T c ^ it is no longer possible to neglect 
the nonlinear contribution which causes m to relax as a function of time according to a power 
law. This general scenario can be modified by the presence of confining walls, which change the 
spectrum of the operator —A in Eq. (|147|) . Thus, if this spectrum has a lower bound above zero 
(as it is the case for Dirichlet-Dirichlet boundary conditions in a film for which the zero mode 
is not allowed) the linear regime eventually dominates even at bulk criticality £ = oo. This is 
a consequence of the critical-point shift in the film geometry. On the other hand, for T < T c ^, 
£~ 2 has to be replaced by — £ _2 /2 (within MFT) and so, for T = T C (L) < T c ^, the spectrum 
of the operator —A — £~ 2 /2 includes and the nonlinear contribution is no longer negligible. 
Therefore at T C {L), m relaxes according to a power law as a function of time. Beyond MFT this 
power law is characterized by the critical exponents of a d — 1-dimensional bulk system. 

Before discussing the nonlinear relaxation in the confined geometry, we summarize briefly 
the results of the corresponding analysis of the semi- infinite geometry, i.e., the effects of a single 
surface [54-56]. In particular we consider the typical relaxation process, realized by applying 
an external field h(x±,t) = h{x±)9{—t) where x± is the normal distance from the single wall. 
In this case the subsequent evolution depends on the resulting initial order parameter profile 
mo(x±) = m(x±,t = 0). Within the field-theoretical approach it is possible to compute the 
scaling function for the evolving order parameter profile, as carried out in Refs. [54] and [56] for 
Model A: 

m(x ± ,t;r) = m Q {t/T Q y^{{xJi G ){t/T Q )- l l z , t/T R , {(m (x' ± )/m )(t/T )' 3 ^}) (148) 

where the temperature enters via the relaxation time Tr (see Eq. ©)• For (mo(Xj_)/mo) 
x (t/To)P/ uz 3> 1 the behavior of the system becomes independent of the initial configuration and 
follows a universal scaling function ^(v, w, oo). From a short-distance expansion it is found that 
^(v — > 0, io, oo) ~ that at criticality the magnetization close to the surface decays 

as t~^ x l vz for t — * oo. For fixed x± and at criticality, after some (non-universal) initial transient 
behavior due to the finite initial magnetization, the order parameter relaxes as m ~ t~^' vz , i.e., 
according to the "bulk" behavior. (Of course one has \&(oo, oo) = VPbulk^))- As time passes 
(i.e., t/To > (x±/^o) z ) the effect of the surface reaches the point x± and the relaxation crosses 
over towards the surface behavior, i.e., m ~ t^^ 1 l vz . (Note that for the ordinary surface uni- 
versality class 0i > (3.) This crossover is nicely displayed in Monte Carlo simulation data [55]. 
Off criticality this picture is modified only by the fact that when m becomes sufficiently small, 
i.e., t 3> Tr, the system enters the linear regime where the decay of m is finally exponential. 
Moreover the influence of the surface penetrates into the bulk only for a finite distance ~ £ from 
the surface. Well inside the non-critical bulk no crossover is expected between surface and bulk 
relaxation. 

Focusing now on the film geometry, the corresponding nonlinear evolution equation for 
m(x, t) is still given by Eq. (|147f) together with the boundary conditions in Eq. (|13|) . In or- 
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Figure 16: Late relaxation of the order parameter m(x±,t) at bulk criticality in a film of 
thickness L as a function of the reduced time variable t = (i/To)(£o/ L) z , for various values of 
the distance x± = x±/L from one wall. Already for i > 0.2, fh(x±,t ) exhibits, independently of 
x±, its asymptotic long-time behavior, characterized by the exponential decay which is indicated 
as a straight line in the figure. 

der to proceed we introduce dimensionless quantities via 

m(x,t) ee M L m(5LL^m = 4 m( " L ' m(L/e ° )Z) (149) 
V 3! Q m 

(using Eq. (JUJ) with x G x [0, 1]. Thus in the absence of the external field Eq. 1)1471) turns 
into 

<%m(x,i) + [-A x + f + m 2 (x,i)]m(x,i) = with m(x s ,i) = 0, (150) 

where f = (L/£) 2 for r > and f = — l/2(L/£) 2 for r < 0. (We recall that, within mean- 
field theory, £(r — > + ) = r ^ 2 whereas £(r — > 0~) = (— 2ro)~ 1 ^ 2 .) Using the scaled variables 
x, t, and f there is no explicit dependence of the profile fh(5t,t;f) on the size L of the layer. 
Static solutions of Eq. (|150|) have been discussed in detail in the literature for various boundary 
conditions (see, e.g., Refs. [57-59] and references therein). For f > r c (r c < determines the 
critical-point shift) the asymptotic solution of Eq. (|150|) . m^x) = lim.f_ too fh(5c,t), is ?fi 00 (x) = 
0, whereas in the low-temperature phase, i.e., f < f c < 0, the order parameter profile is non- 
trivial and its analytic expression can be found in Appendix[DJ The nonlinear partial differential 
equation (|150j) can be solved numerically. Starting from an arbitrary order parameter profile 
rn-o(x) ee m(x, i = 0), with mo(xg) = 0, the profile evolves according to Eq. (|150|) . In analogy 
with the results for the semi-infinite geometry [54-56], we expect a non-universal behavior in 
the early stage of relaxation due to the fact that the order parameter profile we start with fulfills 
the Dirichlet boundary conditions and thus takes on small values near the confining walls, while 
universal properties are observed in the scaling limit of mo(x) being infinitely large. 

We start with discussing the behavior of the system at bulk criticality f = 0. In Fig. ^] 
we show the results of the numerical solution of Eq. (jl5()(l . evolving from an initial profile that 
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is constant in the directions parallel to the confining walls and has a trapezoidal shape in the 
transverse direction, fulfilling the boundary conditions: m(xj_, i = 0) = Kx±/D for < x± < D, 
fh(x±,i= 0) = K for D < x± < 1 - D, and m(x±,i = 0) = K(l - x±)/D for 1 — D < x± < 1. 
The parameters D and K have been suitably chosen to ensure the stability of the numerical 
solution of the equation (typical choices are D ~ 10 -2 and K ~ 10 2 ). Under this assumption 
the problem depends on a single space variable, given by the distance from a wall < x± < 1/2, 
and on the time variable i. The solution for 1/2 < x± < 1 is obtained by taking into account 
the obvious symmetry of the problem, i.e., fh(x±) = m(l — x±), provided the initial profile is 
chosen to share this symmetry. In Fig. ED different curves refer to different distances from the 
wall. These data demonstrate clearly that asymptotically the relaxation is exponential in time, 
independently of x±, according to the behavior 

m(x,t^oo)~e- 7r2 (*/ To ^°/ i ) z . (151) 

This behavior can be explained by analyzing Eq. (|150|) for the special case we are discussing, 
i.e., 

dtm(x±,t ) + [-dl ± + m 2 (x±,t )]m(x ± ,t ) = (152) 

with fh(x± = 0,i) = fh(x± = 1, t) = 0. In the linear regime (i.e., neglecting fh? in Eq. (|152|) ) 
this equation can be solved straightforwardly, leading to 

oo 

fh(x±,i) = ^a n e _7r n *sm(imx±) (153) 

n=l 

where the coefficients a n are determined by the initial profile. Thus for generic values {a n } 
the leading asymptotic decay is indeed exponential (see Eq. (jl51|0 and its dependence on x± 
is given by sin(7rxj_). As already stated at the beginning of this section, the exponential decay 
of the order parameter is intimately related to the fact that in Eq. (|153|) the sum contains no 
zero mode n = and that we are considering the problem at the bulk critical point, located in 
the disordered phase of the confined system. In the following we consider the behavior at early 
times as shown in the logdog plot in Fig. 1171 

For small values of x± (lower curves in Fig. I17|) we note that the order parameter relaxes 
according to a power law ~ t~ Ks until the crossover towards exponential decay takes place. 
On the other hand, for x± = 1/2 (uppermost curve), the power law is different and relaxation 
follows the law ~ t~ Kb with Kb < n s - For intermediate values of x± there is a crossover between 
the two power laws. In order to elucidate this crossover, we compute the effective exponent of 
the relaxation as the logarithmic derivative of the magnetization profile, i.e., 

it(x±,t) = -j — t • 154) 

A powerdaw behavior corresponds to a value of k independent of i, while an exponential decay 
~ e~ pt leads to k = pi. The behavior of n(x±,t) is shown in Fig. 1181 

In accordance with the relaxation behavior in the semi-infinite geometry we expect the 
following qualitative picture of the order parameter evolution in a film. At the beginning of the 
relaxation process the effect of both boundaries on the relaxation behavior starts to propagate 
into the bulk. At a fixed distance from the wall the relaxation process is characterized by the bulk 
exponent Kb until the influence of the closest surface reaches this point and changes the exponent 
into the surface one k s > Kb- We expect that the final crossover towards the exponential decay 
(due to the presence of two confining walls with Dirichlet boundary conditions) will take place 
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Figure 17: Early relaxation of the order parameter fh(x±,i) at bulk criticality in a film of 
thickness L as a function of the reduced time variable t = (i/To)(£o/ L) z , for various values of 
the distance x± = x±/L from one wall. This is a magnification, on a log- log scale, of the interval 
i < 0.2 in Fig. 1161 The two straight lines indicate the power-law behavior in the bulk (uppermost 
line) and close to surfaces, respectively. The crossover from bulk to surface and from surface to 
film behavior (i.e., exponential decay) is evident. For larger values of x± the crossover occurs 
later. 
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Figure 18: Effective exponent n(x±,i) for the order parameter relaxation in a film of thickness 
L (see Eq. (jl54j) ) as a function of i = (i/To)(£o/ 'L) z and for various values of x± = x±/L. The 
two horizontal dashed lines represent the values of the exponent in the bulk (k = 1/2) and close 
to a surface (k = 1), respectively. For i < 10 -4 the non-universal features of initial relaxation 
are evident. The linear increase of n(x±,t > 1CP 1 ) = pi indicates an exponential relaxation 
~ e~ pt with p = it 2 . 
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via the intermediate stage of surface relaxation only if the spatial point under consideration 
can be reached sufficiently in time by the effect of one wall before the effect of the second wall 
arrives. This picture is confirmed by Fig. 1181 which provides a quantitative analysis. The two 
exponents k are readily determined as k s = 1 and kj, = 1/2, corresponding to 0\/vz and f3/vz, 
respectively, as in the case of the semi-infinite geometry. There is a large degree of freedom in 
the definition of times at which crossovers take place. We define the following typical times: 



i c ,b(x±) 
i c ,bs(x±) 



K{t c>b ,x±) = 0A, (155) 
i^(t c ,bs,xj_) = 0.75 , (156) 
K(t c>e ,x ± ) = 1.1 . (157) 



The first one, i c ^, is a measure of the time required to relax from the initial condition. Its 
dependence on x± carries non-universal information about the specific initial profile. In the 
scaling limit of infinite initial magnetization (in the bulk) t cb = 0. The time i c ^ s is a measure 
of the time required to cross over from the 6ulk to the surface behavior. The value chosen 
in its definition is simply half way between Hf, and n s . The corresponding plot is given in 
Fig. (a). It is useful to remark that i c ^ s {x±) is expected to be finite (given that it is 
still related to "surface" behavior) also in the limit L — > oo at fixed x±, i.e., in the semi- 
infinite geometry. As a consequence, based only on scaling arguments, we can deduce the 
behavior of i c fi s {xi_) for small x±. Since k is dimensionless one expects the scaling behavior 
K(xj_,t,L) = (xj_/L,(t/To)(^o/L) z ) (compare Eq. (JSJ)). Thus, according to the definition 
Eq. (17561 . it follows 

^(^) Z = F cM {xjL). (158) 
A finite non-trivial limit for L — > oo for fixed x± exists provided 

F c , bs (y ^ 0) = C bs y z (159) 

with a constant C bs which implies 



so that 



t c ,bs(x±) = C bs x z x , (161) 



with z = 2 within MFT. This is in agreement with Fig. ^3 (b) (apart from very small values of 
x±_ which are numerically still influenced by the initial relaxation) yielding ~ 0.5. 

The time f c>e measures the time required to enter the linear relaxation regime. As can be 
seen from Fig. 1201 t c ,e(%±) attains a nonzero value for x± — ► 0. This means that in the limit 
L — ► oo at fixed x± one has 

t c ,e(x±,L > x±)/T = D e {L/i ) z , (162) 

with D e ~ 0.076. This divergence of the crossover time as function of L is expected because in 
the semi-infinite geometry, i.e., for L — ► oo, the crossover towards an exponential decay never 
takes place. 
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Fi gure 19: (a) Crossover time t c t b s (see Eqs. (|156jl and (|161|0 between bulk-like and surface-like 
behavior of relaxation, (b) Log-log plot of the crossover time i c ^ s . The straight line indicates 
the quadratic law expected to hold for small x± within mean-field theory according to scaling 
arguments (see Eq. (jl61|0 . For i < 10 -2 there are deviations due to the non-universal initial 
relaxation. Beyond mean- field theory t c ^ s {x± — > 0) ~ Xj_. 
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Figure 20: Crossover time t C)e (see Eq. (|157|l ) towards the linear relaxation behavior characterized 
by an exponential decay of fh(x±,i). 

6 Summary 

Within the field-theoretical approach, we have investigated various universal aspects of the 
relaxational critical dynamics [Model A, Eqs. ©-©] i n hhn geometry with Dirichlet boundary 
conditions, corresponding to the so-called ordinary surface universality class. We have obtained 
the following main results. 

(1) (i) We have provided general scaling properties (see Sec. El and Subsec. 14.1)) for the 
two-point response R and correlation functions C in the film geometry (Fig. [2). In partic- 
ular, their Fourier transforms in the directions parallel to the confining walls (with conju- 
gate momentum p) scale as R(p,xi±,x 2 ±,t) = (o^/T^ 1 )(L/^) 1 ~ ri ~ z ^±(p = pL,x±± = 
xi±/L,x 2± _= x 2 ±/L,i =_(t/T±)(£±/L) z ,L = L/0 and C(p, x 1± , x 2± , t) = (6£/T±) 
{L/£^) l ~ n< £±{p, xi±, x 2 ±, t, L), respectively, in terms of the film thickness L, the bulk 
correlation length £(r = (T — T Cj fe)/T Cj 6 — > ± ) = Co = l r l _l/ ) an d the relaxation time 
Tr(t — * ± ) = T |r|~ above and below the bulk critical temperature T c b, with the 
universal ratio Tq /Tq = 3.3(4) in spatial dimension d = 3 (Appendix EJ). The semi- 
infinite limit of is discussed in Eqs. ()59)) - (j(i2)) and (|91j). The explicit forms of the 
associated universal scaling functions and ^+ have been computed within Gaussian 
approximation (denoted by ^). 

(ii) The time evolution of the mean-field scaling function Sfffl [see Eqs. 1)85)1 and (|18Hj) ] 
for p = and T = T c ^ is shown in Fig. |2] for different values of the scaling variables 
Xi± and i. The curves in Fig. |2 provide (apart from an amplitude) the time evolution 
of the order parameter profile across the film at T = T c ^ after a laterally homogeneous 
perturbation confined to the plane x± = x±± and ~ 5(t) has been applied. A qualitative 
feature of interest is the time ii{x\±) at which the inflection points of these profiles reach 
the boundary x 2 ± = 0, leading to a change in the shape of the profiles (see Fig. |3 and 
Appendix [E|) • 

(2) (i) We have discussed the dynamical effects of an external field /i(x = (xii , x±),t) (conju- 
gate to the order parameter) on the fluctuation-induced Casimir force Fim ( x || > r > A ^ {M x > 01) 
acting at the point Xii on the feft and right confining walls. Its scaling behavior is given 
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by F, (r) (x||,T > 0,L,t,{h(x,t)}) = L- d ^ } (x|| = X|,/L, L, i, {(L/( )^[h(x, t)/f) ]}) 
(see also Eq. H10U|0 , Within Gaussian approximation we have provided the expressions 
for the scaling function J^ dy ^ for two specific instances of external field: (W) h(x, t) = 
hw5(x± — X\±)8(t — 1\) and (P) /i(x, t) = hp5(x^ — x 1 ||)(5(a;_|_ — xi±)5(t — ti). In both cases 

^ (dy) is the sum [Eqs. (jl()7fl and of the siatic Casimir force J^( st ) corresponding to 

the ordinary-ordinary surface universality class considered here, and a dynamic term which 
is quadratic in the scaling variables h\y,p [see Eqs. (jl()8j) and (|114|) ] proportional to the 
amplitudes hw,p of the external field. In particular we focused on the force at T c ^ where 
^( st ) = (d— 1)A < [Eq. jnnj]- For both cases W and P we have studied as function of 
time the maximum A^' P of the field-induced contribution [#^ dy ^ — J^^]/h'^ v p . 

(ii) It turns out that is attained at i = ij(x\±) discussed in the previous point 
(i). Figure 0] shows the dependence of on the normal distance x\± at which the 
perturbation has been applied (see Appendix IF.2|) . Figure El illustrates the time depen- 
dence of the normalized part of the field-induced Casimir force (x\±,i) = [J^/ dy ^ — 
^^}/[h 2 w A^(x 1L )} which decays ~ e" 2 ^'. 

(iii) A^ is attained at i = tM(<5x|| , x\±) (Fig. EI) which depends additionally on the scaled 
lateral distance 5xy = (xn — Xx||)/L between the action of the force and the epicenter of 
the perturbation. The maximum of the force on the walls spreads with an asymptotically 
constant radial velocity which decreases with the film thickness (Eq. (|117j) ). The force 
decreases monotonically for increasing lateral distances |<Jxm|. Figures [7| and |H] show the 
dependence of the corresponding maximum A^(5x.pX\±) (attained at i = tM(<5x|| , x\±)) 
on |<5x||| and xu_. The time dependence of the normalized part of the Casimir force 

&l (5x||,£i_L,t) = [J^ (dy) - &( st) ]/[h 2 P A%(5x\\,x 1± )] is reported in Fig. M it decays ~ 

p(d-l) e -27T 2 i_ 

(3) We have computed the universal scaling function of the dynamical two-point cor- 
relation function [see the previous point (1)] within Gaussian approximation. In Fig. ^] 
the time evolution of c ef^ [see Eq. (|122j) ] for p = and T = T c ^ is shown for different 
values of the scaling variables and i. Figure fTTI refers, instead, to the cases in which 
t = l/(p 2 + L 2 ) is finite, i.e., if T > T cf> or p / 0. 

(4) (i) In view of the possibility to probe the spatial structure of correlations by means of 
neutron or X-rays scattering under grazing incidence, we have investigated the behavior of 
the frequency- and momentum-dependent correlation function in planes parallel to the sur- 
face of the film, i.e., of C(p,xi± = x±,X2± = x±,u) = OQ(L/^) 1 ~ ri+z 'rf±(p,x±,x±,u) = 
uT±(L/£y,L) [see Eq. (Mj)]. 

(ii) Near the walls the surface behavior ^+(p, x± — ► 0, x± — > 0, Co, L) = x 2 ^ 1 ^^ Ue £w(P, 
[see Eq. (|63)) ] is recovered. The various asymptotic behaviors of 'tfyy have been discussed 
in Subsec. O [see Eqs. (fB3 |l -([7T |) and (fl3^ |) -(fT3^ |) .(fT37 |) ]. Within Gaussian approximation 
the scaling function is shown and discussed Figs. ^] and EH I n contrast to the semi- 
infinite geometry, in films 'tfyy attains a finite value at the origin of the (^~ 1 , p, cj)-space 
which diverges for L — > oo. 

(iii) The correlation function in the middle of the film is characterized by the scaling 
function ^/(p,(D,Z) = ^+(p, 1/2, 1/2, Q, L). Its asymptotic behaviors are discussed in 
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Subsec. I4~T1 [see Eqs. (jTSJ-dHOJ) and (tTHl - (tTlH|) . (fT13|) ]. Figures [H and show the shape 
of the scaling function ^ along the axes of the p, w)-space. 

(iv) Comparing the scaling functions ^JS' and it turns out that, although their 
shapes are similar, the latter exhibits more rapid algebraic decays for L, p, lo — > oo. 

(5) (i) We have considered the nonlinear relaxation from an initially ordered state by solving 
numerically the evolution equation for the scaled order parameter profile m(x, i) [see 
Eqs. (|149ft and (|150j) ]. which is proportional to the time-dependent mean value of the 
order parameter ((p(x,t)) across the film. In particular we have analyzed the relaxation 
at the bulk critical point T = T c j, from an initial profile that is laterally constant and 
has a symmetric shape in the transverse direction which fulfills the Dirichlet boundary 
conditions. The universal aspects of the relaxation process are independent of the actual 
shape of the initial profile. In view of this symmetry the scaled order parameter profile 
m(x, t) depends only on x± and i. 

(ii) In Fig.E3the late stage of the relaxation of the order parameter is shown as a function 
of t for various values of For all x±, fh(x±,i — ► oo) displays an exponential decay due 
to the critical point shift in the film geometry. On the other hand, during the early stage 
of relaxation, i.e., for t < 0.1, a crossover between an algebraic surface- and bulk-like 
behavior clearly emerges (see Fig. I17j) . Close to the surface one observes fh(x±,i) ~ 
t-h.K»*) (fr./(vz) = 1 within mean- field theory), whereas upon moving inside the film 
m(x ± ,i) ~ (f3/(vz) = 1/2 within mean- field theory). This crossover is clearly 
detected by the time dependence of the effective exponent K(x±,i) = —dlogfh(x±,t)/di 
shown in Fig. 1181 We have also defined (Eq. H156|) ) and studied fFig. I19|) the typical time 
ic,bs(x±) at which the crossover from the 6ulk to the surface behavior takes place for a given 
distance x± from the near wall, and that one for the crossover to the ultimate exponential 
decay i c>e {xA_) (Eq. (fTST!) and Fig. l20j) . 

(iii) Thus the nonlinear relaxation of the order parameter in film geometry is character- 
ized by a cascade of four clearly identifyable, successive decay modes: nonuniversal initial 
relaxation dominated by the initial profile, bulk-like power-law decay ~ t~@/ uz , surface- 
like power-law decay ~ t~^/ uz , and exponential decay ~ exp(— ir 2 t). The initial relaxation 
lasts, in the present case, up to i < 10 -5 and the exponential decay prevails for t> 10 . 
The crossover time between the two power laws depends on the distance of the point of 
observation from the near wall. 
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A Universal dynamic amplitude ratios 

In this appendix we determine the universal amplitude ratio Tq /Tq as introduced in Sec.|2J We 
consider the model described in Sec. 01 in the absence of any confining wall (bulk behavior) and 
for the specific case of a one-component field (i.e., Ising universality class, iV = 1). The linear 
response of this model in the presence of a finite external field h, and thus of a non-vanishing 
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magnetization m, has been studied in Refs. [50] and [51] within a loop expansion up to one loop. 
The results obtained there allow one to compute the relaxation time of the system (a) in the 
high-temperature phase (r > 0, with h = and thus m = 0) and (b) in the low-temperature 
phase (r < and thus m = Mq ^ also for h = 0). In particular the value of the spontaneous 
magnetization Mq for r < and h = can be determined from the equation of state reported 
in Eq. (5.1) of Ref. [50]: 



M (t < 0) 



(163) 



where e = 4— d, (3 = (1 — e/3)/2+0(e 2 ) is the critical exponent of the spontaneous magnetization 
Mo, r = (T — T Ci b)/T Ci b, jx is a momentum scale (see below), and u* is the fixed-point value of the 
renormalized coupling constant (whose actual value is irrelevant for the rest of the computation). 
The true correlation length in the high-temperature phase is given by 



i + + «(«») 



(164) 



with v = (l + e/6)/2 + 0(e 2 ). Accordingly, // can be expressed in terms of the nonuniversal 
amplitude $ as fi = ^[l + e/12 + 0(e 2 )]. 

The linear response function can be computed from Eq. (5.2) in Ref. [50], taking into account 
the different values of Mq in the cases (a) (Mo = 0) and (b) (see Eq. (|163|0 , This equation 
expresses the deviation m\{t) of the magnetization from the constant value Mo in terms of the 
small magnetic field h\(t) that gives rise to it, and thus provides an expression for the linear 
susceptibility x defined through 



mi{t) 



ds x(t — s)hi(s) . 



(165) 



It is straightforward to express the Fourier transform x±(uj) = J + °°dt e lujt x±{t) of x±(t) f° r 
r ^ as 



X+V) =UJ o 



-*— + 1 - ~ 
^0 6 



+ 0(e 2 ) 



and 

where 
and 

with z = 2 + 0(e 2 ) and 
Fr(x) = 



UJ, 



3 2 
lo q = fi / u 2 (-2r) 1 



+ 0(e 2 



d« 1 + V e "e- 2 "/- 



1+1 A -i jK 1 -! 



(166) 

(167) 

(168) 
(169) 

(170) 



which has a branch cut on the real axis for x > 2. The kinetic coefficient in Eqs. (|168j) 
and ()169|) can be expressed in terms of the non- universal amplitudes T + ac or T + exp , and ^q" 
through, c.f., Eq. JHU) or {HSJ. For instance, fi = (? + CXD ) 2 /T + CXD [l + 0(e' 2 )]. 
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As it is the case for the (bulk) correlation length, there are different possible definitions of 
the relaxation time Tr. For example it can be defined as the inverse characteristic frequency 
(see, e.g., Eq. (3.17) in Ref. [2]), i.e., 



Tfi,ac(r) - 



X _1 (w = 0) duj 



(171) 



uj=0 



where the subscript ac stands for "autocorrelation". Indeed using the FDT (see, e.g., Eq. (|34[l) 
one can show that TR iac coincides with the so-called integrated autocorrelation time of the 
magnetization: 



-Rac 



DC 



where Cm(£) is the two-time dynamic (auto)correlation function of thermal equilibrium fluc- 
tuations (around Mo) of the magnetization; it is given via the FDT (see, e.g., Eq. (|32j0 by 
Cju(i) = fu\°°dsx(s). Another possible definition of Tr introduces the "true" relaxation time 
determined by the position of the frequency pole in x(u;) closest to the real axis: 

r 1 (^ = -^ 1 cxp (r)) = 0. (173) 

From Eq. (JT75|) it follows that T R , 

exp governs the long-time exponential decay ~ sxp( t/T ReX p^ 
of the correlation and linear response function away from T c . 
Using the previous expressions one finds that 

Tr _ J KT 1 (1 + b) + 0(e 2 ) , r>0, 

TR ^-\ H)-Hl-le)+0(e>), r<0, (174) 
which yields the universal amplitude ratio 

T °^ ac - 2 VZ (l + L) + 0(e 2 ) . (175) 



rp- 

J 0,ac 



Instead, for Tr ex p one finds 



T„ -i (^r 1 (1 + ^ + 0(e 2 ) , r>0, 

J-R,exp — <\ ti -\-\ fi , fl ln2\,l , rM ,2\ \ 110 ) 



K-)- 1 [l + (i-^) e ]+0(e 2 ), r<0 
so that the corresponding universal amplitude ratio is given by 

Sep = 2 VZ (l + l ^e) + 0(e 2 ) . (177) 

^O.exp V 2 / 

In Ref. [40] the purely relaxational dynamics of the Ising model has been studied both 
analytically and numerically (in three dimensions). In particular the generalizations of Tfj jac 
and T R . 

exp t° the case of a non-vanishing wavevector have been considered for a generic point in 
the (r, Mo)-plane, and their universal scaling functions have been computed (see Eqs. (15), (16), 
and (19) therein). These results provide predictions also for the universal ratios /Tq^. 
and r ~ exp /T ~ ac . Using the notations of Ref. [40], T^ p /T^ ac = 7^ xp (0;oo) = 1 + 0(e 2 ) and 
T o7ex P / T o7ac = ^exp(0;-l) = 1 + (3/8 - In 2/2) e + 0(e 2 ), respectively. It is straightforward to 
check that the results reported in Eqs. (|174|) and (|176j) are in accordance with these predictions. 
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B Static correlation function 



The static (i.e., equal-time) correlation function cffl can be obtained from Eq. (j30j) by integrat- 
ing C (0) (q n ,Lj) over uj (see Eqs. (j2U) and ((23)1): 

C£ ) (p,x l± ,x 2 ±) = J2^ n (x 1± ;L)^ n (x 2 ±;L)- T —- . (178) 

n=l Qri + r ° 

Using the relations (see Eq. (|18j)) 

$ n (xi;L)$ n (x 2 ;L) = - |cos[— (xi - x 2 )] - cos[— (xi + x 2 )]| , (179) 
= P 2 + n 2 n 2 /L 2 , and the identity (see, e.g., §5.4.5-1 in Ref. [60]) 



-^s^ = ri cosh[a(N -!)] _! | for _ 2<z < 2) (18Q) 



n=l 



7r 2 n 2 + a 2 \2a sinha 2a 2 J 



(outside the range \z\ < 2 this expression has to be extended periodically in z with period 2), 
one finds 



(o), , 1 cosh[a(|xi_L - x 2 ±\ - L)\ - cosh[a(xi_L + x 2 ± - L)\ 

C y J (p,x 1± ,x 2 ±) = . h j \ ( 181 ) 

2a smh(aL) 



where a 2 = p 2 + £ 2 (see footnote [SJ) . By applying addition formulae for hyperbolic functions, 
Eq. (|181|) can be expressed as 

r (°)( n r r \ - T sin HaL(x</L)} sinh[aL(l - x>/L)} 

<^st {P^i±,x 2 ±) - L . , (182) 

aL sinh(aL) 

where x^ = max{xu_, x 2 ±} = (x%± + x 2 ±_ + \x\±_ —x 2 ±_\)/2 and x^ = min{xij_, x 2 ±} = (xu_ + 
x 2 i_ — \x%± — X2_i_|)/2, in agreement with Eq. (B4) in Ref. [44] (and with Eq. (4.1) of Ref. [6] in 
the case cq = +oo). This result can also be found directly from Eq. (|95|) by using the FDT (see 
Eq. ©): C£\p,x 1± ,x 2± ) =C<®(p,x 1± ,x 2 ±,t = 0) = R^(p,x 1± ,x 2X ,u = 0). 

C Response function and its asymptotic time dependence 

With a view of the scaling function ^ for the response function (see Eqs. (|82[) and (|83|l) we 
consider the sum 

oo 1 

£ $n(xi;L)<S> n (x 2 ; Lje-^"/^ 1 = -V(x x /L, x 2 /L, nt/L 2 ) (183) 

n=l 

with the eigenfunctions <I> n (x;L) defined in Eq. (|18jl. Using Eq. (|179(> . Eq. (|183|) reduces to the 
evaluation of expressions of the form 

OO 11 + 00 

^cos(7rn^)e-" 2n2T = -^ + - e ^-^nV (lg4) 

n=l n=— oo 
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which can be expressed, according to standard definitions (see, e.g., §16.27.3 in Ref. [61]), in 
terms of Jacobi's theta function 



+00 



M*,q)= E v n e2mz 



(185) 



so that with t = £lt/L 2 (see Eq. l|81jl) and Xi = Xi/L 



Xl-X 2 _ n 2£ Xl+X 2 _ 7r 2 f 

#3^ ,e n ) — i9 3 (7r ,e 71 ) 



(186) 



By using Poisson's resummation formula [i.e., ^2'^L_ 00 e 2mnx = YlnL-00 ~ n )\ we obtain 



i} 3 (irx,e 



-TV t\ 



•x) a /f 



(187) 



so that 



^f(x 1 ,x 2 ,t) 



'Ant 



e -(xi-x 2 ) 2 /(4t) _ e -(x 1+ x 2 ) 2 /(4i) 



(188) 



+ 



E [^ 



e v 2 



e v 2 



n=l 



This expression is useful for discussing the semi-infinite limit of the response function (see 
Eq. ©). 

The long-time limit t = Qt/L 2 — > 00 of Eq. (|183j) follows from 



tf 3 (7rx, e _7r *) = 1 + 2e~ 7V £ cos(2vrx) + 0(e 



(189) 



so that 



*(xi,x 2 ,i -> 00) = e _7r *[cos7r(xi - x 2 ) - cosvr(xi + z 2 )] + 0(e ) (190) 



as expected also from Eq. (|183|) because in that limit the sum in Eq. (|183|) is dominated by its 
first term so that 



V(xi,x 2 ,t ->• 00) = L$i(xi; L)$x(x 2 ; L)e 77 + 0(e 



(191) 



D Analytic expression for the mean-field order parameter pro- 
file 

The rescaled mean- field order parameter profile rn-oo (x) is the stationary solution of Eq. Q15U|) : 

[-Ax + r + m^ (x)]m 00 (x) = (192) 

with the Dirichlet boundary conditions m^x^) = 0. Thus ^(x) is an extremum of the 
functional (03), which can be rewritten as 



H[m] = L d - 4 m§^ / d d x 



i(Vxm(x)) 2 + ^rm(x) 2 + ^m(x) 4 



(193) 
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by taking into account Eqs. (|5i|) and (|149j) . the definitions x = x/L, r = {L/^) 2 = L 2 for r > 0, 
f = -l/2(L/£) 2 = -Z 2 /2 for r < 0, and V = R d-1 x [0, 1] as well as m(x B ) = 0. We recall 
that, within mean-field theory, £(r — » + ) = r whereas £(r — > 0~) = (— 2ro) -1 / 2 (Eq. (JHJ)); 
in the main text we use also the abbreviation L/^ = L. Here and in the following £ means for 
r < and for r > 0. In view of the translational symmetry in all directions parallel to the 
confining walls, finding the equilibrium profile reduces to solving a one-dimensional differential 
equation for ip(x±) = m(x) where x = (R\\,x±): 

+ fip(x ± ) + V 3 (x ± ) = , 
V(0) = ^(1)=0, 

with the symmetry ip(x±) = ip{\ — x±) so that for regular solutions ip'(l/2) = 0. Accordingly a 
set of equivalent boundary conditions for Eq. I)194j) is ip(0) = and ip'(l/2) = 0. Equation (|194j) 
has always the trivial solution ip(x±) = 0. Note that Eq. (|194ft describes the one-dimensional 
closed motion of a particle with coordinate ^ in a potential V(tp) = -fip 2 /2-ip 4 /4 function 
of "time" x±. Using this analogy one can show that the non-trivial solution vp(x±) of Eq. (|194j) 
is bounded by the non-trivial solution ip of Eq. H194|) with free boundaries resembling the bulk 
solution, 

* - { °^ £ I < o ; ^ 

in the sense that |"0(icj_)| < 4> 6 . In view of this, in the following we consider f < 0. Equa- 
tion (|194j) can be integrated once by using the boundary condition ^(0) = 0: 

?P' 2 (x±) = ^(0) + f^ 2 (x ± ) + ^\x ± ) . (196) 

The solution of Eq. (|196|) (together with the boundary condition ^(0) = 0) depends on ^'(O) 
which has to be determined as a function of f in order to fulfill the second boundary condition 
ip'(l/2) = 0. It is convenient to introduce the scale transformation ij)(x±) = V2k( a(C,x±), with 
er(0) = and cr'(0) = 1. In terms of <j(w = Ca?j_) Eq. Q196|) turns into 

<r /2 = (1-<7 2 )(1-A;V) , (197) 

where 

-f = ( 2 {k 2 + 1) . (198) 

The two parameters k > and ( > just introduced have to be determined as a function of f 
in such a way that <r'(£/2) = 0, corresponding to the condition tp'(l/2) = \/2&C 2<j/ (C/2) = 0, 
and that Eq. (|198|) is satisfied. Equation (|197jl is solved by the Jacobian elliptic integral of the 
first kind with modulus k (assuming that ip'(x±) = ^/2k(, 2 a'((x±) > for < x± < 1/2), 

<r(w) ^ 

ds - = w , (199) 

o ^(l-s 2 )(l-k 2 s 2 ) 

and thus (see, e.g., 8.144.1 in Ref. [62]) 

a{w) = sn(w; k) . (200) 

6 For f < 0, ip is the position of the maximum of V(tp) for ip > 0. A particle starting from the point ip = at 
the "time" x± = with tp'(x±) > cannot perform a closed motion (and then return to the starting point at the 
"time" ii = 1) if its initial kinetic energy is large enough to overcome the potential barrier V(ip), which allows 
for ip > ip. As a consequence, for a closed motion one has \ip \ < ip. 
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The parameters k and £ are related via Eq. (|198|) . They are fixed by imposing the boundary 
condition o'(Q/2) = 0. Thus Eq. (j!97j) renders two possible values: cr((/2) = 1 and <r(C/2) = 
1/k 2 . On the other hand, Eq. (fT93j) implies th at V>(l/2) = y/f = V2k(a((/2) < $ = CVk^+T 
(where we used Eq. flSE^), so that <r(C/2) < ^(/c 2 + l)/(2/c 2 ), i.e., (a) a((/2) = 1 for < k < 1 
and (b) a{Q/2) = 1/k 2 for > 1. We first consider case (a). Using a(Q/2) = 1 in Eq. Q199|) one 
finds the following relation between £ and k (see, e.g., 8.111.2 and 8.112.1 in Ref. [62]): 

K(k) = | , (201) 

where i"T(&;) is the complete elliptic integral of the first kind. This allows one to replace the 
variable ( in Eq. lfT9*5J) : 

-f =[2K{k)] 2 {k 2 + 1) . (202) 

This is an implicit equation k = k(f) for the modulus in terms of the physical variable f = 
-l/2(L/£) 2 (as defined for f < 0). Thus the solution of Eqs. (jUg) and (fTMj) is given by 

tp(x ± ) = 2\[2kK(k) sn(2K(k)x ± ; k) . (203) 

Since K(k) is a monotonicly increasing function with i^(0) = vr/2, one has [2K(k)] 2 (k 2 + 1) > it 2 
and therefore there is a non-trivial solution ip(x±) ^ only for 

r < -vr 2 . (204) 

This determines the well-known mean-field critical point shift in the film geometry with Dirichlet 
boundary conditions on both sides [57,58,63,64] (y = 1/2): 



-^(f) 1/l/ - (205) 

The discussion of case (b) proceeds accordingly. Using the property sn(ku;l/k) = ksn(u;k) 
(see, e.g., formula 106.01 in Ref. [65]) one finds that the solution for case (b) is the same as for 
case (a) provided the modulus k is replaced by 1/k. 

These results can also be obtained from the findings in Ref. [59], in which the mean- field 
order parameter profile is analytically determined for the case of (+,+) and (+, — ) boundary 
conditions (see Ref. [59] for details). Alternatively, Eq. (|194|) can be integrated so that 

4>' 2 (x ± ) = fV 2 (x ± ) + ^\x ± ) - f^ 2 (l/2) - ^ 4 (l/2) , (206) 

using ip'(l/2) = which is valid both for Dirchlet-Dirichlet and (+, +) boundary conditions. 
The function u(x±) defined as 

satisfies Eq. (|206[) . provided that A is chosen according to 

^ = -fV 2 (l/2) - 1^(1/2) , (208) 

which is equivalent to — rti 2 (l/2) — l/2u (1/2) = A 2 /2. Therefore both ip(x±) and u(x±) 
solve Eq. ()206|) . If ip satisfies (+,+) boundary conditions, then u satisfies Dirichlet-Dirichlet 
boundary conditions and vice versa. Thus it is possible to take advantage of the result reported 
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in Ref. [59] for the former case in order to solve the latter we are presently interested in. (Note 
that the normalizations used here are different from those used in Ref. [59], see Eq. (|194j) here 
and Eq. (A3) therein.) In Ref. [59] two different parameterizations are provided for the solution 
denoted as m++: one for tL 2 > —it 2 in Eq. (A13) therein, and the other for tL 2 < —it 2 in 
Eq. (A15) therein. One can see that in the former case the value of A 2 is negative, i.e., the 
corresponding solution for Dirichlet-Dirichlet boundary conditions is not real, and one is left 
only with the trivial solution which is identically zero. Instead, in the latter case, taking into 
account the different normalizations, 

= *-rmk) m(iK \ k)z . k) , (209) 

where k is determined according to Eq. (|2U2|) . Thus A 2 = QAk 2 [K{k)] 4 and the corresponding 
function Up^^z) = A/tp +j+ (z) is identical to the solution ip(z) given in Eq. (|2U3[) . 

From Eqs. (|203j) and (|202|) it is possible to recover the result for the mean-field profile in the 
semi-infinite geometry by considering the limit of large L for fixed £_ and x± which corresponds 
to |r| 3> 1, x± <C 1. For |f | 1 one has k — > 1~ and thus one can use the approximation (see, 
e.g., formula 17.3.26 in Ref. [61]) 

K(k) = -\ In (^f-) (1 + 0(1 " *)) • (210) 

Given that sn(ii, 1) = tanh u (see, e.g., formula 127.02 in Ref. [65]) one easily finds at leading 
order 

ij)(x±) ~ ip tanh ^j-^-xj^J for x± <C 1, |f| > 1 . (211) 

Using Eqs. (|149j) and (|195|) and the fact that within mean- field theory Co /Co = V% one can 
express the profile as 

m(x,r) = m %tanh (l^f) — > vaoT l/2 (l - 2e~ x ^-) , (212) 

where here £ = ^(r < 0) = This expression agrees with the well-known result for the semi- 
infinite geometry [66,67]. In Fig. |^ the order parameter profile ip(x±)/il> (normalized to the 
corresponding bulk value t/j, see Eq. (|195|) ) is shown for some values of f. The comparison with 
the profile in the semi-infinite geometry (Eq. (|211|) ) is also shown. 

From Fig. one can infer that for f — > — oo the order parameter profile in the middle of 
the film x±_ = 1/2 approaches rapidly the bulk value. Indeed, defining 5ip = ip(x± = 1/2) — ip 
one finds _ 

- S 4 = l-^, (213) 

where k is determined by Eq. (|202j) . Using Eq. ()210|) one finds that for \f \ 3> 1 (i.e., in the limit 
of large film thickness L at a fixed temperature) 

Jjt = Ae-^ 2 (l + 0(e-^ 2 )) , (214) 

if) 

i.e., in the middle of the film the deviation of the order parameter profile from the bulk value 
due to the distant confining walls decays ~ exp[— L/(2£)]. 
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Figure 21: Order parameter profile ip(x±) (normalized to the corresponding bulk value ijj, see 
Eq. (Q2SJ)) across the film, for some values of f = -l/2(L/£) 2 : -10 (a), -15 (b), -20 (c), -40 
(d), —80 (e). The order parameter vanishes for f > — it 2 ~ —9.87 (Eq. 1)2041) ). The dashed lines 
for x± < 1/2 represent the corresponding profiles (normalized to the bulk value tp) in the semi- 
infinite geometry (given by Eq. (|211|0 in which the order parameter vanishes only for r > 0. In 
the cases (d) and (e) these provide good approximations to the actual order parameter profile 
for x± < 1/2; for (e) the differences are barely visible. 

E Relaxation of the response function 

In the main text tj = Ftj(x~i±) has been defined as the reduced time at which the inflection 
point of the scaling function ^(xi±,X2±,i) for the mean-field response function (as a function 
of X2±, see Fig. |2j) close to x± = disappears. The position of this point x^ = Xj(x\±,i) is 
determined by the condition 

^°' 2) (x l± ,4 ± ,i) = 0, (215) 

where, here and in the following, we use the notation ^>^ n ' m > (xi±, x 2 x, t ) = d£ i± d™ 2± i S>(xi±,X2±,i) 
Accordingly, ii is given by the time at which the inflection point reaches the surface at x± = 0, 
i.e., it is implicitly determined by the condition Xj(xi±,tj) = 0. This provides an equation for 
the function F tl (see Eq. (|9~2|l ): 

X I (x 1± ,F tl (x 1± )) = . (216) 

The function shown in Fig. has been determined by solving numerically Eqs. (|215j) and (|216|) . 
Here we determine analytically the behavior of Ftj(xi±) for x±± — > and x~i± — > 1. From 
Eq. (fT86|) it follows that 

^(xi±,X2±,i) = ^(x2±,x 1± ,i), $(xi±,X2±,t) = -*(*i±, -x 2 ±,i) , (217) 

and 

V(x 1± = 0,x 2± ,i) = *(x 1± = l,xz±,t) = . (218) 

Accordingly ^ ^{x 1± = 0,x 2± ,i) = ^ ' m \ Xll _ = l,x 2 ±,i) = and ^°' 2m Hx 1± ,x 2 ± = 
0,t) = 0. (Eq. H217|) implies that is an odd function of X2± if analytically continued to 
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negative values of X2±_ by using Eq. (|186[> .) Therefore Eq. (|215() is always trivially satisfied for 
X2± = 0, whatever the values of x±±_ and i are. The function Ft { is given by the nontrivial 
solution of 

* {0 ' 2) (x 1± ,x 2± 0,F tl (x 1± )) = . (219) 
This solution can be found by considering the series expansion of f(°' 2 ) for %2±. — * 0, i.e., 

y [0 ' 2) (x 1± ,x 2 ± -> 0,F tl (x 1± )) = ¥°> 3 \x 1± ,0,F tl (x 1± )) x 2 ± + 0(x 2 2± ) . (220) 

Therefore ti = F tl (xij_) is determined by the condition 

^ (0 ' 3) (x 1± ,0,t 7 ) = . (221) 

Let us consider the limiting case x\± — ► 0, i.e., the behavior of -Ft 7 (y — > 0). As discussed in 
Subsec. l4~2*l (see Eq. (|94[)). in this limit one expects — ► 0. Using Eqs. (|186|) and (|187|) . one 
realizes that in the limit i / — * the terms of the sums with n = are the leading contributions 
to ^(°' 3 )(xi_L,0,t/). (The terms with n ^ give rise to exponentially small corrections.) Thus 
Eq. flZZD) reduces to 

(xf ± - 6ti)Si_L = , (222) 
that, apart from the expected trivial solution x~ij_ = 0, is solved by 

-2 

*/ = -M^LL - 0) = . (223) 

The full numerical solution shown in Fig. [31 is in accordance with this analytic result. 

We now consider the case x%± — > 1, i.e., F tl (l). Expressing Eq. ()221|) as before, one realizes 
that the leading contributions to vp( ' 3 ) (xi± — > 1,0, ij) stem from the terms in the sums with 
n = 0,l (Eqs. H186jl and (|187[) ). Indeed Fig. shows that ij is significantly smaller than 1 
also for x\±_ = 1, allowing one to neglect terms that are suppressed by a factor ~ exp(— 
compared to the leading one. This leads to 

x 1± (x 2 1± - 6tj) - (2 - x 1± ) [(2 - x 1± ) 2 - Mi) exp ( Xl± .~ 1 j = , (224) 

which has the expected trivial solution x\± = 1, whereas the nontrivial one for xi± — ► 1 is given 
by 

*/-*J + -£ = 0, (225) 

yielding 

h = F tl (x 1± = 1) ~ 1 _ -L ~ 0.09176. (226) 

This provides a very accurate approximation of the actual value F tl (l). Note that the second 
solution of the quadratic equation (|225|) can be discarded because it is inconsistent with the 
assumption ij <C 1, under which Eq. (|225j) has been derived. The actual value of ^(1) can be 
computed from Eq. (|221|) after having discarded the trivial solution xi± = 1. This can be done 
considering the expansion of ^(°' 3 ) around x±± = 1, retaining only the leading term. This yields 
the implicit equation 

* (1 ' 3) (l,0,t 7 ) = . (227) 

Here we do not report the corresponding expression that can be easily worked out. The corre- 
sponding solution can be numerically determined and is given by 

F tj (l) ~ 0.0917918. (228) 
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F Computation of the Casimir force 



F.l General expression 

According to Eq. (|103|) it is possible to compute the force exerted on the confining walls in 
terms of the stress-tensor Tj_^. In order to determine the effect of a time-dependent external 
field on this force we first note that according to Eq. (J12[) within Gaussian approximation, due 
to {(p)o = and ([y] n [<^] m )o = for m > n, one has 

(<p(xi,ii)<p(x 2 ,£ 2 ))ft = M x i>*iM x 2,*2))o 

+ JdV'dt'dV"dt" h(x! ,t')h(x!' ,t") (229) 

x n 2 (<p(*', £'M*i, ti)}o(0(x", t"Mx2, t 2 ))o 

(higher order terms in h are zero because n = 2) where has been introduced after Eq. 
Using Eq. (|229jl one finds 



(T_L±(x,t))hUedv = \ dx 1± d X2± ((p{x l7 t)ip(x 2 ,t)) h \ xi=X2=x€dv 

I ,■ ^ ( 23 °) 

(T±±)o\dv + 2 



where the response function R is defined in Eq. (|12|). The previous equation provides the 
expression for the force density i*}( r )(x||) acting on the left (right) plate (depending on which 
part of the two boundaries dV is used in Eq. 1)230(1 ) for a general external field. The previous 
expression can be expressed in terms of the corresponding dimensionless scaling variables (see 
Eqs. (@3J), (53), JB5J), JSE), and jHEJ): 

^ (dy)(0) (L,t,{^}) = ^( st )( )(0,X)+i[ / -^^dxLLdti fcfoWl) 

8LJ (2vr) ^ (231) 

x e-^ 2 H d i2X ^ (0) (p, x 1± , x 2 ± , t - h , L)\ S2±=0 

where the scaling variable h is defined by 

h(p,x 1L ,h) = C^^L/^f^L-^hip/^^L^n^o/L)-^ (232) 

with 

h(p,x ± ,t) = J d^^n /i(x||,a; ± ,f)e ip - x il . (233) 

h is the analogue of h introduced after Eq. (|105() . Note that the former carries, compared to 
the latter, an extra factor L~ stemming from the Fourier transform in the d — 1 parallel 
spatial directions. Moreover, within the Gaussian model (go = 0) it is not possible to define, 
in accordance with Eq. (|99j) . the scale factor f)o for the external field h, whose engineering 
dimension £ ( rf+2 )/ 2 follows from the Gaussian action. 

Note that, according to Eq. (|230|) . the dependence of the Casimir force on the applied field 
is quadratic and therefore the effect of a sum of fields hi + h>2 is not equal to the sum of the 
separate effects of each field. 
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F.2 Specific cases 



In Subsec. 14.31 we consider two different instances of externally applied fields: a perturbation 
h(x,t) = h\y5(x± — xi±)5(t — t\) which is spatially constant in the plane x± = x\± parallel 
to the confining walls and a perturbation h(x,t) = hp5(x\\ — Xi\\)5(x± — xi±_)5(t — t±) that is 
localized at a point x = (xjn , X\±) within the film. In the following we present the details of the 
computation of some relevant quantities determining the response to these external fields. 



F.2.1 Planar perturbation 

Here we discuss the asymptotic behaviors of the amplitude A^{xn_) defined in Sec. 14.31 (see 
Eq. which determines the Casimir force maximum (d — 1)A + h^yA^ (x\±). In terms of 

the notation introduced after Eq. (j215j) . is given by 



2 L 



(234) 



Using Eqs. (|186|) and (|187|) one finds an expression of vlK ' 1 ) which is suited to discuss some 
asymptotic behaviors at short times: 



*(°' 1 )(x 1± ,0,t) 



E 



n 



xi± \ 
2 J 



exp 



i / x 1± y 



"t " 



(235) 



Accordingly, 
<Br(o.i)(x 1± -.l,0,f 



(l-xi±) 

+oo 
X E 



1 (n-1/2) 2 

2 t 





" (n-1/2) 2 " 


exp 


i 



(236) 



+ 0((l-x 1± ) 2 ) 



and 



x l± 



0,0,i) = ^-g^ {exp[-z 2 ± /(4t)] +0(e- 1 ' 1 )] . (237) 

(Note that the approximation provided by this expression deteriorates upon increasing i.) In 
Appendix FBI we found that, at leading orders, ii(x~x± — ► 1) = F tl (l) ~ 0.0918 and ii(xx± — > 
0) = zf ± /6. From Eqs. (|234l and (I236J) one has 



with 

and where 



A%(x 1± ^l) = ai (l-x 1± ) 2 
ai = £i (*/(!)) ^ 17.535 



(4?ri ) d - x 2^? 



E 



1 (n-1/2) 2 1 

2 t 



exp 



(n-1/2) 5 



(238) 
(239) 

(240) 



this function will be useful also for the discussion of the localized perturbation below. In the 
opposite limit xx± — > Eqs. (|234|) and (|237f) yield 



Aj(xx± - 0) = /dM ± 



(241) 
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with 



JCi ~ 0.427888 



where we have introduced the constant 



K,,] = — 
2 



1 [2(d + 2)] 



d+2 



-(d+2) 



(47T) 



(242) 



(243) 



Moreover, the asymptotic expressions for ^I^ ' 1 ) in Eqs. (|236f) and (|237|) provide the corresponding 
ones for the function (see Eq. pi2j) ) because 



For xi± — > 1, Eq. (|236j) gives 



£l(fr(l)) ' 

whereas, for xu_ — ► 0, Eq. (|237|) and ij(xi_L — ► 0) = x 2 ^/6 render 



- I e 

s 



-3/a 



S = F/tj(x 1± -► 0) . 



(244) 



(245) 



(246) 



These asymptotic behaviors are indicated in Fig. [5J Equation (|115|) allows one to determine the 
long-time behavior of <^]± (xi±,i) for fixed xi±: 



^(x 1± ,i^ oo) ~ e- 2 ^ f [l + 0(e- 37r2f )] 



(247) 



which clearly displays the expected exponential decay for t 3> l/(37r 2 ). 



F.2.2 Localized perturbation 

In Subsec. POl we have introduced the amplitude A£(5xm 
response to a point-like field (see Eq. I|12U|) ): 



xii — x.i\\,xi±) associated with the 



1 e -(5x||) 2 /[2tM(5S||,Six)] f 

2 [4^*^(^X11, xij.)] 



d-1 



* (0,1) (a;i±, *Af ((5xi| , x 1± ))\ S2±=0 



(248) 



Here we discuss its asymptotic behavior for |<5x||| — > 0. The typical time iw(5xn, jcijJ when the 
effect of the perturbation attains its maximum at a point with \Sx.»\ <C 1 is expected to be of 
the same order as that in the case of a planar perturbation at the same xi±_, which is given by 
the function ti(xi±) shown in Fig. |31 Accordingly, fjvf(5xii — * 0,x\±) < 0.1 and therefore, in 
Eq. (|248p we can use the expressions given in Eqs. (|236|) and ()237j) for ^C 3 ' 1 ) in order to obtain 
the asymptotic behaviors for x\± — ► 1 and xu_ — > 0, respectively. We expect that, as it is the 
case for in Fig. 0J they provide good approximations of the actual dependence. Taking into 



account that ijv/(<5x|| 
one finds 



and in particular 



0,x 



i± 



0) 



+ (5x,,) 2 ]/[2(d + 2)] (Eq. (dHJ)) and using Eq. (12371) 



a£ (ten 



0,s 



1_L 



0)=ACd 



[x 2 ± + (5X||) 2 ] d + 2 ' 



AX(<5x,| = 0,xi± 



0) = /Q/^ +1) 



(249) 



(250) 
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For d = 3 one has /C3 ~ 0.169773 whereas /C4 ~ 0.148406. For x\±_ — ► 1 we consider only the 
case 5x|| = 0. From Eq. (|236[) . one finds 

A£(5x|| = 0, xll -»■ 1) = a rf (l - Xll f , (251) 

where (see Eq. (|2l 



a d = Ajfaf (<**|| = 0, s 1± -► 1)) . (252) 

Taking into account that [see the discussion after Eq. l|118Jl ] tjv/(^ x || = 0,xu_ — ► 1) = (4 — 
y/U)/10 ~ 0.0683375 for d = 3 and t M {&x.\\ = 0,x 1± -»■ 1) = (9->/§7)/20 ~ 0.0604236 for d = 4, 
one has 03 = £3(^(^x11 = 0,xu_ — > 1)) ~ 17.927 and 04 = C^ijai^w = ^>^1-L — * -0) — 19.416. 

Let us discuss the behavior of A^(6xn,Xx±) for large |5xn | and < x\± < 1 fixed. According 
to Eq. (|116|) one has, in leading order, i^^XM — ► oo,xij_) = |<5xn |/(27r) — (d — l)/(47r 2 ) + 
0(l/|(5x|||) and therefore one can use in Eq. (|248|) the asymptotic behavior of <I'( ' 1 ) reported in 
Eq. (HE}, finding 



e 



-27r|5xii| 



A£(5x|| -> 00, x 1± ) = 27T 2 ^^_ sin^vrxn.) . (253) 



For |5x|||,xij_ — ► one obtains from Eq. |2 

A%(6x.u,xi± 



A£(Jx|| =0,xi ± ) [l + (,5x||)2/x 2 J d+2 



(254) 



which indeed provides a very good approximation of the actual curves already for \Sxn | , < 0.5 
(see Fig.©. 



G Expansion of the in-plane correlation function 

In Subsec. [4.4l we discuss the mean-field expression for the correlation function for points located 
within a plane parallel to the confining walls, i.e., C^°\p, x±, x±, lj), with the corresponding 
scaling function given in Eq. ()127l) . This expression depends actually on two variables: a = aL 
and = xj_/L (0 < x± < 1). In order to discuss the behavior of Eq. (|127|) for \d\ <C 1 (so that 
\a(l — 2x±)\ <C 1), we write 



cosh a — cosha(l — 2x±) 
a sinh a 



J2^n(x ± )d 2n (255) 



n=0 



where V n is a polynomial of degree 2n + 2, has the symmetry V n (x±) = V n {\ — x±), and 
T'n(O) = 0- From the series expansion of the l.h.s. of Eq. ()255() one easily finds that 

V (x) = 2x(l - x) , (256) 
Vi(x) = -^x 2 (l-x) 2 , (257) 

V 2 {x) = —x 2 (l-x) 2 (l + 2x- 2x 2 ) , (258) 
Vs(x) = -^x 2 {\ - x ) 2 (l + 2x-^-- 3x 3 + ^ . (259) 
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Taking into account that a? = p 2 + I? — iu> (see Eq. (|97|l). one finds that 
cosh a — cosha(l — 2x±) 



Im- 



a sinh a 

oo 2fc+l<n / s (260) 



n=l k=0 

and thus 

&°\p,x ± ,x ± ,u,L) 



En^) E (- 1 ) fc+1 ( 2 , n + i)(p 2 + I2 ) ? 



o<j 2fc+l<n ✓ s 

2E^(^/i) E (-i) fe+1 L n , 1 )(p 2 +^^ l2,+,! 11 



(261) 



n=l fc=0 
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